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FINAL TECHNICAL REPORT - NAGW-68 
Studies of Planetary Upper Atmospheres Through Occultatlons 

The goal of our research Is to understand the structure, 
composition, dynamics and energy balance of planetary upper 
atmospheres, through proper Interpretation of stellar 
occultatlon data. Stellar occultatlons are sensitive probes of 
the structure of planetry upper atmosphere, but it has not 
proven easy to extract definitive results from the 
observations. The spiky structure of the light curves has been 
variously interpreted as scintillation due to Isotropic 
turbulence and as evidence for layered atmospheric structure. 
Correct identification of the physical mechanisms at work will 
provide essential information for understanding the dynamics 
and energy balance in the upper atmospheres of Uranus, Neptune, 
and Jupiter. Our main objects, under this grant, have been: 

1. To investigate the wave-optical problem of modelling 
strong scintillation for arbitrary turbulent 
atmospheres. 

2. To determine whether turbulence dominates the 
observations . 

3. To relate these results to the physical processes in 
the observed upper atmospheres. 

In the appendices, we have included manuscripts which have been 
submitted for publication, and which contain the detailed 
results of our investigations. Here, we shall summarize 
briefly the principal results of our studies. 

We began our work with a rigorous investigation of the 
effects of turbulence on stellar occultatlon light curves, to 
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learn how the observations could be used to learn the relative 
Importance of turbulence and other dynamical processes In the 
regions probed. To study these questions, we developed a 
computer code capable of generating occultatlon light curves, 
using wave optics, by specifying the refractive properties of 
the model atmosphere. We confirmed the analytic weak 

scintillation theory of Hubbard ^ al. (1978), but found ^ 

significant qualitative differences under conditions of strong 

scintillation. Our key finding was that Kolmogorov turbulence, 

whether Isotropic or anisotropic, can account for neither the 

details of the spiky structure in the light curves nor the 

large-ampl t tude temperature variations found by numerical 

Inversion of the data. This was the first definitive test of 

Young's (1976) contention that the observations are dominated 

by the effects of scintillations. 

We applied our newly developed analysis techniques to the 
Uranus occultatlon of 15 August 1980. Observations separated 
by many Fresnel scales along the limb of the planet show a high 
degree of correlation, suggesting that the atmosphere Is 
strongly layered. We confirmed that the derived variations in 
the thermal profile are not artifacts of Kolmogorov turbulence. 

Next, we investigated the energy balance in the Uranus 
upper atmospliere, and considered meclmnisms that might maintain 
the large obsiM'ved tempera tv; re variations against a host of 
disruptive processes. We found evidence that the mean 
ttunperature of the upper amtosphere has changed significantly 
over a period of four yeai.:, and investigated the associated 
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energetics. We concluded that molecular and eddy diffusion, 
and atmospheric dynamics, are as Important as radiation to the 
upper atmospheric heat balance. 

Finally, we investigated the properties of light curves 
produced by a wide range of realistic atmospheric models. 
Including fields of Inertia-gravity waves in the presence of 
turbulence, and determined the effects of severe ray crossing 
on the numerical inversion method of obtaining temperature 
profiles. We concluded that quasi-global features of 
atmospheric structure are accurately determined by numerical 
inversion; horizontally inhomogeneous structures (including 
turbulence) with length scale L << (2irRH)^/2, where R ■ 
planetary radius and H » scale height, are filtered out and 
have little effect on the derived profiles. 

We are now in a position to extend our work to detailed 
modelling of light curves under the specific conditions 
appropriate to individual occultation events, with the goal of 
comparing the upper atmospheric structure of the outer planets. 
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1. The Upper Atmosphere of Uranus : A Critical Test of 

Isotropic Turbulence Models , by French, R.G., J.L. Elliot, 
B. Sicardy , P. Nicholson, and K, Matthews, Icarus (1982), 
in press. (Appendix B), 

2. The Thermal Structure and Energy Balance of the Uranian 
Upper Atmosphere , by French, R.G., J.L. Elliot, E.W. 
Dunham, D. Allen, J.H. Elias, J.A. Frogel, and W. Liller. 
Icarus (1982), submitted and refereed; currently under 
revision for final s’ibnission. (Appendix C). 

3. Strong Turbulence and Atmospheric Waves in Stellar 
OccultatioDS , by French, R.G. and R.V.E. Lovelace, Icarur. 
(1982), submitted. (Appendix D). 

4. Turbulence Models and the Atmosphere of Uranus , by French, 

R.G. and J.L. Elliot, Bull. Am. Astron. Soc. 704. 

Abstract. (Appendix E). 

5. Structure of the Upper Atmo sphere of Ura lus , by French, 
R.G., April, 198T lAU/RAS Colloquium on Uranus and the 
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Observations of the 15 August 1980 Uranus occultation of KM 12. obtained from Cerro Tololo 
InterAmerican Observatory, European Southern Observatory, and Cerro Las Campanas Observa- 
tory, are used to compare the atmospheric structure at points separated by ~140 km along the 
planetary hmb. The results reveal striking, but by no means perfect, correlation of the light curves, 
ruling out isotropic turbulence as the cause of the light curve spikes. The atmosphere is strongly 
layered, and any acceptable turbulence model must accommodate the axial ratios of 260 which are 
observed. The mean temperature of the atmosphere is 150 — I5*K for the region near number 
density 10'^ cm~’. Derived temperature variations of vertical scale **130 km and amplitude ^5*K 
are in agreement for all stations, and correlated spikes correspond to low-amplitude temperature 
variations with a vertical scale of several kilometers. 


I. INTRODUCTION 

Observations of planetary upper atmo- 
spheres are of interest because they pro- 
vide information about the radiative and 
dynamical processes at work in these rare- 
fied regions, and permit tests of models of 
thermal structure, atmospheric wave prop- 
nation, and turbulence. With the advent of 
high-speed photometry, stellar occultation 
observations have been made in recent 


years of the atmospheres of Jupiter, Nep- 
tune, Mars, and Uranus [see Elliot (1979) 
for an extensive review and references]. 
The occultation light curves almost invari- 
ably contain sharp intensity variations, or 
spikes, indicating that the atmospheres are 
not isothermal, but there has been no agree- 
ment on the nature of the atmospheric re- 
firactivity variations responsible for the 
spikes. 

Two approaches to the problem have 
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been adopted. The hrst has been to assume 
that the atmosphere can be {q>proximately 
characterized by homogeneous layers of 
large horizontal extent, that geometric op- 
tics is applicable, and that ray crossing is 
not severe. Under these circumstances, a 
unique mapping exists between the ob- 
served light curve and the atmospheric tem- 
perature prohle. The alternative strategy 
has been to regard spikes as manifestations 
of turbulence-induced scintillation. Ac- 
cording to this view, in its most extreme 
form, the assumption of spherical symme- 
try is not even approximately satisfied, and 
the detailed temperature structure deduced 
from numerical inversion of a light curve 
bears little relation to the true atmospheric 
basic state. 

In the past, there has been no definitive 
test of these opposing views. Careful stud- 
ies of the observations of the Jupiter occul- 
tation of 3 Sco reached fundamentally dif- 
ferent conclusions about the importance of 
turbulence and the structure of the Jovian 
atmosphere (Elliot and Veverka, 1976: 
Hunten and Veverka, 1976; and Jokipii and 
Hubbard, 1977). Combes et al. (1975) 
pointed out that the divergences among the 
inferred temperature profiles are the results 
of different analysis methods, but there re- 
mains a clear difference of opinion about 
the importance of turbulent scintillation in 
stellar occultation data. From multiple ob- 
servations of the Mars occultation of c 
Gem, French and Elliot (1979) concluded 
that the numerical inversion method was 
valid for Mars, and that the temperature 
profiles provided evidence for semiglobal 
waves, similar to those observed by Viking 
1 and Viking 2 (Seiff and Kirk, 1977) and 
predicted by Zurek (1976). However, the 
Texas- Arizona Occultation Group (1977) 
and Hubbard (1979) found the observations 
equally consistent with isotropic turbulence 
theory. As a final example, large-amplitude 
fluctuations were observed in the Mariner 5 
Venus occultation. These data provide a 
demonstration of the small effect on tem- 
perature profiles of these fluctuations. 


based on the separate frequency and ampli- 
tude inversions of Mariner S data which 
yielded nearly identical profiles (Fjeldbo et 
al., 1971). 

It is important to resolve the dispute. F 
the occultation temperature profiles are re- 
liable. they provide information about the 
dynamical state of the upper atmospheres 
of Mars, Neptune, Jupiter, and Uranus. On 
the other hand, if the effects of turbulence 
are dominant, it may be possible to com- 
pare the turbulence fields of these planets. 
"A sensitive test which can in principle de- 
termine which hypothesis is correct would 
be to observe the same occultation from 
two points on Earth for which the ray paths 
are separated by - 100 km along the limb. A 
significant correlation would favor the lay- 
ering hypothesis, whereas a lack of correla- 
tion would favor the turbulence model" 
(Jokipii and Hubbard, 1977). In this paper 
we report the results of such a test. On 15 
August 1980. Uranus occulted star No. 12 
on the list of Klemola and Marsden (1977). 
High-quality photometric observations at 
2.2 p.m wavelength were made from several 
stations, and the results indicate striking, 
but by no means perfect, correlation of 
structure in the Uranus atmosphere for dis- 
tances of ~140 km along the planetary 
limb. We conclude that the atmosphere is 
strongly layered, and that the observed 
spikes are not manifestations of isotropic 
turbulence, or of any reasonable form of 
anisotropic turbulence. 

II. OBSERVATIONS 

In order to test turbulence models, we 
will compare the occultation light curves 
obtained from Cerro Tololo Inter American 
Observatory (CTIO), European Southern , 
Observatory (ESO), and Cerro Las Campa- 
oas Observatory (CLCO). Several charac- 
teristics of the observations are given in Ta- 
ble I. Here we restrict our attention to 
information provided by a direct compari- 
son of the data. Further analysis and inter- 
pretation, and detailed information about 
each observation, can be found iiHjie re- 
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TABLE * 

OCCULTATION ORSERVATIUNS 


L :iser- 
vstory 

Telescope 

aperture 

(m) 

Filter 

References 

cno 

4.0 

K 

French ei al. (1982) 

ESO 

36 

K 

' 'cardy ei at. 1 1982) 

CLCO 

-.5 

K 

Niclioisoii ff 'll. 
(1982c) 


ports of each g;.-un. cited in Table I. All of 
the obscp ’-ifO'’* • , made in the K band 

fXfl i 7 .itr. L,\ ^ 0.4 am), where the 
irom strongly depressed. 

' i.e ih''iC stations san',*.*” J regions sepa- 
rated by up to “140 km along the limb of 
r anus. Figure 1 shows the geometry of the 
event. Immersion occurred at a latitude on 
Uranus of approximately -26°S and probed 
the dark limb of the planet; emersion oc- 
curred near latitude + 18®N and probed the 
illuminated limb. All stations obtained high- 
quality photometric data of the immersion 


event, and 300 sec of the corresponding 
! light curves are shown in Fig. 2. A strong 
correlation of several of the large intensity 
fluctuations is immediately apparent, as 
will be shown quantitatively in the next sec- 
tion. 

Each group experiemed difficulties dur- 
ing the emersion phase of the occultation: 
at CnO, clouds caused variations in signal 
level; at ESO, problems with guiding 
caused the star to drift in and out of the 
^rture; and at CLCO, recording was not 
begun until mid-emersion. In spite of these 
problems, convincing correlation of several 
large spikes in the emersion light curves is 
also evident (Fig. 3). 

For this analysis, the CLCO data were 
digitized from a copy of the stripchart rec- 
ord of ;he observations. Nicholson et al. 
(1982a) independently digitized the ori^nal 
record at a higher time resolution for their 
analysis. Slight differences in the digitiza- 
tions are responsible for some of the differ- 
ences between the CLCO results reported 
here and those given by Nicholson et al. 
The inherent difficulty in the digitization 


IMMERSION EMERSION 



Fio. I. OcculuUon feoinctry. The immersion nod emersioa inckt of the ttnr relative to the plane- 
tary limb are shown for each station. Immersioa occurred near latitude - 26* uid emersion occurred 
near latitude 18*. The symbols on each track show the appamt stellar location as seen simuluneously 
from each station at a time prior to immersion and at a time after emersion. 




ORIGINAL PAGE IS 
OF POOR QUALITY 
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Fio. 2. Immereion occulution Ught curvet. Time it in lecondt after the following reference times on 
IJ Aug I9M (UT): CTIO. 22* 36- 5«0I; tSO. 22* J7- 02t40; CLCO, 22* J7- 03«65. Thete timet were 
chosen to that the tight curvet would be in appnnimate alignment at the attrometric half-light time. 
The sgiproximate altitude scale in the atmosphere relative to the half-hgin level it shown. The cross- 
conelatioo results shown in Fig. 4 were computed over the interval labeled CCF. Unit stellar intensity 
Icveb are shown on the left and zero imensity levels on the right. 
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procedure is in maintaining a constant zero 
level for each section of the chart record, 
which in our analysis is uncertain by about 
1% of the stellar intensity. Direct digital re- 
cording of the data onto magnetic tape was 
not possible at CLCO, but would have elim- 
inated this source of error. 

111. AI4ALYSIS AND DISCUSSION 

The extent of similarity in atmospheric 
structure along the limb can be determined 
by examining correlations in the raw data 
and by comparing temperature profiles ob- 
tained by numerical inversion. The most di- 
rect test is cross-correlation of the light 
curves themselves. We begin our analysis 
by computing the intrinsic rms spike ampli- 
tude, a quantity which can be directly re- 
lated to turbulence theory, and estimating 
the effects of noise. We then discuss occul- 
tation timing and proper alignment of the 
light Ciirves for the subsequent cross-corre- 
Uuion analysis, and we compare the corre- 
lation results with predictions of turbulence 


theory. On the basis cf this evidence of ex- 
tensive layering in tlte Uranus atmosphere, 
we then adopt the assumptions of the nu- 
merical inversion method and compare the 
temperature structure derived from the 
three separate immersion observations. Fi- 
nally, we discuss the relation of spikes to 
the derived temperature profiles. 

A. Comparison of Light Curves 

Amplitude of light curve variations. A 
quantitative measure of the strength of the 
Donisothermal character of a light curve is 
provided by the mean-square amplitude of 
the spikes. In the context of turbulence the- 
ory, this is known as the scintillation 
power, and can be dir^tly reUded to the 
amplitude of the putative turbulent refrac- 
tivity fluctuations. To determine the intrin- 
sic strength of the light curve variations, it 
is necessary to detennine the contribution 
of noise to the observed spike amplitude. 

To estimate the importance of noise, we 
computed the quantity c, the >ms noise in 
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Fig. 3. Emersion occultation ti^ht curves. Time is in seconds alter the following reference tunes on 
16 August 1980 (UT): CTIO. 0^ 15" 38!2l: ESO. 0* 15" 3lt40; CLCO. 0^ 15" 29!45. Thew times were 
chosen so that the light curves would be in approsimaic alignment at the astrometric half-lighi time. 
The approximate altitude scale in the atmosphere relative to the half>lighl level is shown. The cross* 
correlation results shown in Fig. 4 were computed over the interval labeled CCF. Zero rtellar intensity 
levels arc shown on the kit and unit intensitv levels on the right. 


the light curve tor a I -sec integration, ex- 
pressed as a fraction of n«, the counts per 
second for the unoccuJted star. If rt( ij) is the 
mean counting rate in counts per second for 
the ^h integration bin of duration It, N is 
the total number of integration bins used in 
the calculation, and n is the mean counting 
rate in counts per second for all N integra- 
tion bins, then 


(A/)''* ( I A 

■ ‘"'1 • 


Using a total of 50 sec of data (at tin i rso- 
lution Af “ 0!2) prior to immersion atW af- 
ter emersion, when the only variations in 
the light curve would be due to noise, and 
computing the average c for a running mean 


of 4 sec duration, we obtained the results 
given in Table II. The principal uncertainty 
in the ESO emer ion results is in determin- 
ing n*, because the star was drifting in and 
out o( the aperture. (This was true to some 
extent during immersion as well, since otT- 
set guiding was not possible until after im- 
mersion had occurred.) llie quantity ( vj 
is a dimensionless measure of the rms 
noise of an occultation (French ei at., 
1978), where v. is the velocity of the event 
perpendicular to the limb of the planet and 
M is the mean scale height. We computed t „ 
from event astrometry and assumed that H 
■■ 70 km. The noise is lowest for the CTIO 
observations, whkh us^ the laigest-aper- 
ture telescr^. 
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TABLE li 


Light Cuitvt Fluctvationi 



Evi.it 

V* 

(km scc''> 

t 

(V./WV 


cno 

Inunenion 

S.lll 

0.007 

0.002 


ESO 

ImmcnkM 

6.792 

O.OI? 

0.004 

O.IS 

CLCO 

Imincnion 

6.788 

0.0I3 

0.004 

O.IS 

cno 

EmenkNi 

6.888 

0.008 

0.003 

0.28 

ESO* 

Emenion 

6.874 

(0.U4) 

(O.0I3) 

(0.45) 

CLCO 

Emersion 

6.8S8 

O.OI I 

0.003 

0.32 


* Perpendiculu' velociiict computed for an oblate planet with an ellipticity <rf’ 0.024 (Elliot cr cl., 1981). 

* A Kate height of 70 km was assumed, cotresponding to a tempenturc u( I3I*K. 

' Results for the ESO emersion dau are affected by guiding problems during the gradual reappearance of the 
occulted star. 


Having determined the contribution of 
noise to the light curve fluctuations, we 
computed the rms amplitude of 

light curve fluctuations due to the atmo* 
sphere of Uranus, from 

where A<t>< is the ith observ ed residual of the 
light curve from a running mean of 4 sec of 
data, normalized by the unocculted stellar 
intensity. For the regions of Figs. 2 and 3 
marked **CCF," we obtained the results in 
Table 11. In all cases, the observed fluctua- 
tions are far larger than the noise level. 
Photon noise is not a significant source of 
cnor compared to known systematic ef- 
fects, such as guiding errors and variable 
atmospheric transparency. The combined 
eflects of all noise sources air substantially 
smaller than was the cas'* for the p Sco oc- 
cultation by Jupiter (Jokipii and Hubbard. 
1977) or the c Gem occultation by Mars 
(Texas- Arizona Occultation Group, 1977; 
French and Elliot. 19/9). (For example, the 
terrestrial scintillation power was approxi- 
mately flve times the Martian scintillation 
power in the McDonald Observatory < Gem 
data analyzed by the Texas-Arizona Oc- 
cultation Group.) The current observations 
are the least noisy of any which have been 


used to study atmospheric correlations 
along pla. *t'ary limbs. There is good agree- 
ment among stations for the intrinsic spike 
amplitude during immersion » 

0. 16). and although the emersion results are 
less certain, the rms spike amplitude is 
about twice that during immersion. 

Occultation timing and light curve align- 
ment. In order to compare the light curves 
in defail, it was necessary to determine ref- 
erence times at which the occulted star was 
at a common sky-plane distance from the 
center of Uranus for all stations, both for 
immersion and for emersion. It has become 
conventional to use the “half-light" time as 
a reference: this is the time at which the 
best-fitting model isothermal light curve 
falls to half-intensity. We computed tin, fhe 
isothermal fit half-light time, for ail the im- 
mersion data aiKl for the dlO emersion 
light curve. The results of these fits are 
given in Table ill. Because of variable sig- 
nal level and insufflcieni data for the three 
emersion light curves (see Fig. 3), we have 
not included a comparison of all of the 
emersion half-light times. Nicholson et al. 
(1982a) discuss the emersion results in 
more detail and relate them to the determi- 
nation of the oblateness Uranus. 

An independent set a( reference times 
was comiHited from the event astrometry. 
The sky-plane coordinates of the center of 
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TABLE in 


OCCULTATION TIMINQ RESULTS 


Nation 

Event 

Isothermal fit 
f,5 (UT) 

Astrometric* 
fn (UT) 

(sec) 

cno 

1 

IS Aug 1980 

22* 38- 09107 

(22* 38" 09!07) 

(0.00) 

ESO 

1 

IS Aug 1980 

22 38 09.02 

22 38 13.49 

-4.47 

CLCO 

1 

IS Aug 1980 

22 38 13.34' 

22 38 14.79 

-I.4S 

cno 

E 

16 Aug 1^ 

0* 16" 26142 

((F 16- 26142) 

(0.00) 

ESO 

E 

16 Aug t980 

— 

O 16 19.60 

— 

CLCO 

E 

16 Aug 1980 

4 

O 16 17.65 

— 


* The astrometric half-light time is the time at which the occulted star was at the sky-plane distance from the 
center of Uranus corresponding to the CTIO isothermal 6t half-li^t time. The effects refraction were assumed 
to be the same for all three stations. 

* The difference between the isothermal fit and astrometric half-light times. 

' Nicholson ei at. (1982a) obtained an isothermal fit half-light time of 22** 38* 14^ from their analysis of the 
CLCO immersion data, with a corresponding Ari>] of 0.0 sec. 

* Nicholson et al. (1982a) obtained a half-light time of O’* 16" 07*J. based on an isothermal fit of the CLCO 
emersion light curve in Fig. 3. The discrepancy between the CTIO and CLCO emersion isothermal half-light 
times is probably caused by the combined effects of incomplete dain for CLCO and variable transmission for 

cno. 


the planet were determined to high accu- 
racy from a simultaneous solution for the 
orbits of the nine known Uranian rings, us- 
ing all of the available rir,g occultation data 
(Elliot et al., 19811. We assumed that the 
planetary center corresponds to the center 
of mass of the ring-planet system. With 
this information and the known relative 
sky-plane velocity of the planet and star, 
we determined an “astrometric half-light 
time." This is the time at which the oc- 
culted star was at the SKy-plane distance 
from the center of Uranus corresponding to 
the CTIO isothermal ht half-light time. 

If the atmosphere were horizontally ho- 
mogeneous and the light curves were noise- 
free, then the astrometric and isothermal fit 
half-light times would be identical. In fact, 
there are diflercnces of up to several sec- 
onds, as seen in Table 111. To estimate the 
effects of the known noise level on the iso- 
thermal half-light times, we superimposed 
preimmersion and postemersion data on an 
isothermal light curve and compared the 
best-fit half-light times to the noise-free val- 
ues. We found that the known noise intro- 
duced an uncertainty of about ±0!8, larger 


than expected on the basis of photon noise 
alone (French et al., 1978). but less than the 
largest Af|/7 in Table 111. Nicholson et al. 
(1982a) obtained a CLCO immersion half- 
light time corresponding to A(|n = 0.0 sec 
from their analysis of the CLCO data. The 
laige discrepancy between the astrometric 
and isothermal fit half-light times for ESO is 
probably a result of low-frequency modula- 
tion of the light curve due to guiding errors 
during immersion, since offset guiding was 
not emirioyed. Comparison of the light 
curves in Fig. 2 reveals a loss of signal in 
the ESO data relative to the other two ob- 
servations over the interval of 100 to 130 
sec, consistent with this interpretation. 
Henceforth, unless otherwise noted, w'e 
shall use the term “half-light time" to refer 
to the astrometric half-light time. We em- 
phasize that these times are computed rela- 
tive to the CTIO isothermal fit times. 

In the past, isothermal fit half-light times 
have been used for astrometric determina- 
tion of planetary oblateness for Uranus (El- 
liot et al., 1980, 1981), Mars (French and 
Taylor, 1981), Neptune (Kovalevsky and 
Link, 1969), and Jupiter (Hubbard and Van 
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Flandern, 1972; Lecb'acheux et aL, 1973; 
Taylor, 1974). This is a legitimate proce- 
dure as long as it is rect^ized that tlw half- 
light altitude may not correspond to an iso- 
baric or a constant density surface. 
However, the present results suggest that, 
occasion, isothermal half-light times 
may be more uncertain than is suggested by 
the level of noise in the data. The mutu^ 
consistency of half-light times can be 
checked only with multiple observations of 
the same event. The ESO immersion iso- 
thermal fit and astrometric half-light times 
(Table III) differ by more than 4 sec, corre- 
sponding to an error in the deduced radius 
of v^A /]/2 - 30 km. Because of the orienta- 
tion of the Uranus pole, occultation obser- 
vations are restricted to the equatorial re- 
gions. The total difference in planetary 
radius at the equator and at 30° latitude is 
only -- 150 km, so that this newly identified 
potential source of error should be consid- 
ered when estimating the uncertainties of 
the oblateness determination and the im- 
plied rotation period. 

Cross-correlation of light curve resid- 
uals. One measure of the similarity of atmo- 
spheric structure '• .ing the limb is the ex- 
tent of agreement in the detailed features of 
the light curves among the three observa- 
tions. Although the relationship between 
light curve spikes and atmospheric struc- 
ture is indirect, it is useful to quantify in a 
direct manner the degree ot simiiarity uf the 
observations themselves before proceeding 
to numerical inversion of the light curves. 

We computed the cross-correlation func- 
tion (CCF) of the light curve residuals for 
aU of the immersion and emersion observa- 
tions. Our method is similar to that used by 
French and Elliot ( 1979) in a comparison of 
occultation observations of Mars. Residual 
curves at 0!2 resolution were formed by 
subtracting the running mean of 4 sec of 
data from each point; the results were es- 
sentially unchanged when 8-sec running 
means were employed. The data interval 
used was similar for both immersion and 
emersion: from 10 sec above to 20 sec be- 


low the half-light level of each curve (see 
Figs. 2 and 3). The corresponding ran^ in 
normalized stellar flux is Ixtween 0.63 and 
0.31. The residual curves were aligned at 
the astrometric half-light times, and the 
CCF was computed for lags between ±5 
sec in 0!2 increments. Each curve was 
scaled in time by the appropriate prior to 
performing the cross-correlation. 

The results are shown in Fig. 4. For both 
immersion and emersion, there is a con- 
vincing peak near zero lag, and the shape of 
the CCF is similar for both sets of three 
curves. The maximum observed correla- 
tions, the corresponding lags, and the sepa- 
ration along the limb between each pair of 
stations are given in Table IV. Also in- 
cluded are the predicted time oflTsets be- 
tween the astrometric half-light times. For 
a uniformly layered atmosphere on a 
spherical planet, the maximum correlation 
would occur at zero lag. The observed lags 


CTIO-CLCC 

immersion C'lO-ESO 


z eXO-ESe 

O i 



Fio. 4. Immersion and emersion cross-correlation of 
liaht curve residuals. Sec (ext for details of the compu- 
utioo. 
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TABLE IV 

Cross-Comelation Results 


Station 

Event 

d- 

Cl..,* 

.Astro- 

metric 

lir 

iscci 

Las'* 

(seel 

CLCO- 

ESO 

1 

t: 

0.7? 

I..T0 

-0.3 

cno - 

ESO 

1 

ii: 

0.67 

- 4.47 

0.0 

CTIO- 

CLCO 

1 

144 

O.W 

-5. ’’2 

+ 0.1 

CLCO- 


ESO 

E 

33 

10.841 

-I.V5 

-O.I 

CTIO - 

ESO 

E 

115 

I0.7.T) 

6.8: 

-0.1 

CTIO - 

CLCO 

E 

147 

0.10 

8.77 

00 


* </ is the distance atong the Uranus limb separating the 
regions sampled by the two stations. 

* < «w is the obsersed maximum correlation. Results invoK - 
ing the ESO emefsioa data are affected by guiding ptoMcms 
during the gradual reappearance of the occulted star 

' Astrometric Ai is the predicted time offset between the 
astrometric halMighl times of the two statioas. in the sense 
first station - second station. 

The lag corresponding to r™, is: observed time offset of 
the second station relative to the lirsi sution for maximum 
correlation - astrometric At. 


arc pleasingly small. As a coniinnation of 
the smallness of any systematic timing er- 
rors among the three stations, Nicholson et 
al. (1982b) obtained mean discrepancies 
(observed - astrometric) in the ring occul- 
tation times of 0.2 sec (ESO - CLCO) and 
:s0. 1 sec (CTIO - CLCO). The full width 
at half-ma)dmum (FWHM) of the CCF is 
~0!8, or about 5.4 km along the light curve. 
The characteristic FWHM of the CTIO im- 
mersion autocorrelation function (ACF) is 
~0?6, corresponding to 4. 1 km. On the ba- 
sis of models of ring profiles, an upper limit 
to the subteni of the occulted star is esti- 
mated to be 0.18 milliarcsec, or 2.4 km at 
the distance of Uranus. The Fresnel scale, 
(X,D)‘'-. is 2.45 km, where X is the observed 
wavelength and D is the distance to the 
planet. The observed FWHM of both the 
ACF and the CCFs exceed both of these 
length scales, but it is still possible that the 
observed light curve is significantly filtered 
by convclution with the stellar image and 
by the diffraction resolution limit. 

The maximum correlation as a function 
of station separation is shown in Fig. 5. For 
immersion, the correlation remains high for 


separations of over 100 km, decreasii^ 
slowly with distance. FtM* emersion, the 
correlation is even higher, due perhaps to 
the larger intrinsic spike amplitude (Table 
II), but the trend with increasing distance is 
unclear because of the difficulty at ESO of 
acquiring the star in the aperture during 
emersion. 

The rotation of Uranus during the occul- 
tation and the singular viewing geometry 
during the event cause the sampled suboc- 
cultation points in the Uranian atmosphere 
to be separated by somewhat more than the 
distance along the limb jMTOjected into the 
sky plane. For a rotation period of 16 hr. 
the immersion CTIO — CLCO separation at 
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Fig. 5. Cross-corrclAtion of light curve residuals as a 
function of separation distance along Die limb of the 
planet. For immersion, the correlation remains high 
for separations of over 100 km, decreasing slowly with 
distance. For emersion, the correlation is very hi^. 
but the trend with increasing station separation is un- 
clear because of guiding problems during ESO cmer 
Sion. The correlation extends much funber than pre- 
dicted on the basis of isotropic turbulence theory 
(Jokipii and HuMiard, 1977) and provides evidence 
that the atmosphere is strongly stratified. 
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Uranus is increased by 16 km and the emer- 
sion separation by 25 km over the values 
given in Table IV and shown in Fig. 5. 

As a test of the degree of correlation 
throughout the occultation, we determiMd 
the peak cross-correlation immersion 
light curve residuals as a function of posi- 
tion along the light curve. The results are 
displayed in Fig. 6, where we plotted the 
peak cross-correlation as a function of the 
mid-point, relative to the half-light time, of 
the 30-sec interval of data used to compute 
the CCF. The mid-time was stepped in S- 
sec increments, and the maximum correla- 
tion within I sec of the zero lag position 
was recorded. For the interval shown, the 
peak correlation is rather uniform over a 
broad section of the immersion light 
curves. (Simiiar calculations were not pos- 
sible for the emersion data because of the 
problems previously cited.) The results 
shown in Figs. 4 and 5 correspond to a mid- 
time of +5 sec. 

Although the correlation is strong during 


the interval used in the computation, the 
detaUed agreeiiwnt armuig the light curves 
deteriwates de^r into the occultation. 
The immersion data in Fig. 2 have common 
features well after 200 sec, but it is unclear 
how many of them are chance alignments. 
Deep in the occultation, a slight difference 
in depth in the atmosphere, Az. can cause a 
wide separation in time, Ar, between fea- 
tures on the light curve. At low stellar flux 
(A 0. 

A/ “ Az/v>. (3) 

Fm- example, consider a thin spike-pro- 
ducing layer that differs in altitude by only 
2 km over the ISO-km horizontal distance 
separating the repons sampled by two sta- 
tions. At an altitude in the atmosphere of 
200 km below the half-light level, corre- 
sponding to A O.OS, the time separation 
along the light curves is Ar » S sec. Since 
the typical spike separation on a single light 
curve is 10 sec in this region of the light 
curve (Fig. 2c), such slight altitude varia- 



*CCF 


Fig. b. Peak cross<orr«latioo of immcnion Ufhi curve residuals as a function of position akHtg tbc 
li^t curve. The abKissa indkaies the tune, relative to the half-liflit time, trf the mid-point of the 30- 
sec interval of dau used to compute the cross-correlation. 
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lions would make it impossible to correlate 
the spikes in different light curves unambig- 
uously. 

Figures 2 and 3 have been labeled with an 
altitude scale, relative to the half-light 
level, computed assuming an isothermal at- 
mosphere of scale height 70 km and a per- 
pendicular veloc.iV of 6.8 km sec~L (The 
altitude scale is only approximate, since the 
atmosphere is not isothermal.) Wave-opti- 
cal calculations of strongly focused spikes 
(French and Lovelace, 1982) reveal that an 
individual spike can be transformed into a 
cluster of spikes as the overfocused condi- 
tion is reached. We tentatively identify the 
CTIO spikes at 265 and 275 sec in Fig. 2c as 
candidates for this diffraction phenomenon. 

A key diagnostic of the importance of 
turbulence is the degree of correlation of 
spikes in the light curse as a function of 
separation distance. For weak isotropic tur- 
bulence. the CCF FWHM would be roughly 
the larger of the Fresnel scale or the pro- 
jected stellar diameter, and there would be 
no significant correlation for stations sepa- 
rated by distances much larger than the 
FWHM. For the 15 August 1980 event, the 
isotropic turbulence theory (Jokipii and 
Hubbard. 1977) would predict the CCF 
shown in Fig. 5. The observations could 
hardly be in greater disagreement. The cor- 
relation is still strong, even for separation 
distances of 60 vacuum Fresnel lengths. 
Hubbard ei t:l. (1978) have developed a 
wave-optical theory for weak anisotropic 
turbulence, which can be formally extended 
to perfectly stratified layers parallel to the 
limb. However, the present results imply 
an aspect ratio of well over 60 for the turbu- 
lent “pancakes.'* Under these circum- 
stances. the scaling arguments from which 
the Kolmogorov turbulence spectrum is de- 
rived are questionable. As Jokipii and Hub- 
bard (1977) point out: “It is possible to con- 
struct extremely anisotropic turbulence 
models for which the correlation distance 
parallel to the limb is larger than a few kilo- 
meters. but these arc effectively layering 
models." Model realizations of anisotropic 


Kolmogorov turbulence, valid even for 
strong scintillations, show very little resem- 
blance to the observations (French and 
Lovelace. 1^2), and we conclude that the 
effects of Kolmogorov turbulence do not 
dominate the structure of the Uranus upper 
atmosphere. 

The situation may well be different on 
other planets. Woo et at. (1980) found evi- 
dence for anisotropic Kolmogorov turbu- 
lence in the Venus atmosphere on the basis 
of radio occultation measurements, for ex- 
ample. However, based on the Pioneer Ve- 
nus probe measurements, Seiff et al. (1^) 
found that this altitude region has small- 
scale thermal structure (see their Fig. 21). 
Thus it remains possible that layers might 
be important in interpreting the scintilla- 
tions for Venus as well as for Uranus. 

.As a final cautionary note, we mention 
that statistical theories for fluctuations in 
light curves are based upon the starting 
premise of stationary stat'.Uics. A station- 
ary ensemble is obtained only if a large 
number of occultations at identical levels in 
the atmosphere are compared. 

B. Comparison of Temperature Structure 

Mean temperatures and temperature 
profiles. Estimates of the mean temperature 
of the region of the atmosphere sampled 
during the occultation can be obtained both 
by numerical inversion and by modeling the 
light curves. Although the large intrinsic 
spike amplitude guarantees that the atmo- 
sphere is not perfectly isothermal, an iso- 
thermal model light curve fit to the observa- 
tions does provide a measure of the mean 
scale height, which governs the overall rate 
of diminution of stellar intensity. Wc per- 
fonned isothermal fits on all of the immer- 
sion light curves and on the CTIO emersion 
data. The mean temperatures wc obtained 
are shown in Table V. A mean immersion 
temperature of I46"K was found. The CTIO 
emersion temperature was ISb^K. The 
characteristic uncertainty in these values, 
as for the inversion results that follow, is 
±I5'K. 
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TABLE V 


Mean Atmoshiemc Temperatvees 


Sutk» 

Event 

f (botbomal 

m 

(TO 

f (ittvenkm)* 
(*K) 

cno 

1 

131 

136 

ESO 

I 

143 

132 

CLCO 

I 

143 

140 

Avenge 


146 

149 

cno 

E 

136 



*Tbe RM»n inverskm temperatures vc averaged 
over an altitude interval of ~I60 km centered at mim* 
ber density ~3 x lO'* cm'* and pressure -6 x io** 
mbar. 

Using the results of the isothermal fits to 
define the full and zero stellar intensity lev- 
els. we obtained immersion temperature 
profiles by numerical inversion using the 
method of French et al. ( 1978). The method 
assumes that the atmosphere is stratified 
and that ray crossing is not severe. The 
adopted atmospheric parameters used in 
the calculations are listed in Table VI. and 
the profiles are shown in Fig. 7. The profiles 
are aligned in altitude at the level corre- 
sponding to the astrometric half-light time, 
given in Table 111. which for CTIO is at 
number density 1.2 x 10'^ cm~^ and pres- 
sure 3.2 X 10~^ mbar. Because of the insta- 
bility of the initial condition of the inversion 
method, the profiles shown extend only to 
about one scale height above the half-light 
level, and the temperature gradients must 
be regarded as uncertain. The interval 
along the light curve used in the inversion 
can be estimated by comparing the altitude 
scales in Figs. 2 and 7. Similar profiles were 
obtained for the CLCO data by Nicholson 

TABLE VI 


Adopted Atmosphemc Paeametees 


Me*n moleculEr weight 

2.2 gme^' 

Accdcration of gravity 

818 cm sec * 

Atmospheric refractivity (STP) 

1.26 X I0-* 

Tcmpcntura/icEle height 

2.I6*K km-' 



Fk). 7. Temperature pn^les obtamed by nummkal 
iavmkHi of the i m me wto n 1^ curves. The praAies 
are aligned in altitude at the level conespondu^ to the 
astrometnc halM^t time, given in TM)le 111. viitich 
for CnO is at number density 1.3 x |0*^ cm'* and 
pressure 3.2 x |Q-> mbar. 

et al. (1982a). who used the inversion al- 
gorithm of French et al. (1978), and for the 
ESO data by Sicardy et al. (1982), who 
used the inversion technique described by 
Vapillon et al. (1973). We note that Nich- 
olson et al. (1982a) obtained CLCO im- 
mersion temperatures that are as much as 
1S*K lower than our results at the bottom of 
the temperature profiles. These differences 
appear to be attributable to the ~ 1% uncer- 
tainty in the base levels of the indepen- 
dently digitized data. 

If the light curves were noise*free and the 
atmosphere perfectly layered, then the de- 
duced temperature profiles would be identi- 
cal. In fact, althouj^ they have many simi- 
lar features, there are significant dif- 
ferences. The overall temperature gradi- 
ents differ, as do the altitudes of several 
common features. Since the depth in the 
atmos{rfiere is inferred by integrating the 
light curve in time (see French et al. (1978) 
for details], an error of a few percent in the 
assumed zero intensity level can produce 
significant distortions in the temperature 
profile. 

Fortunately, the hi^ degree of correla- 
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tion in the spikes allows us to compensate 
for these distortions. Spikes present in all 
three light curves, displaced in time by the 
differences in astrometric half-light times, 
provide markers which reveal the extent of 
the altitude distortion. In order to compare 
the detailed structure of the temperature 
profiles, we removed the best-fitting linear 
temperature gradient from each profile over 
the altitude range corresponding to the time 
interval from 7 sec prior to astrometric half- 
light to 85 sec after half-light (correspond- 
ing to an altitude range from 25 km above 
half-light to 135 km below the half-light alti- 
tude for the CTIO observations). This is the 
same altitude interval used by Dunham et 
al. (1980) to display the March 1977 Uranus 
atmosphere results. Tlie resulting tempera- 
ture difference profiles ("AT profiles") arc 
shown in Fig. 8 and the mean temperature 
from each fit is given in Table V. The aver- 
age immersion temperature obtained by nu- 
merical inversion is 149°K. Note carefully 



AT (K) 


Fin. 8. Deviations of temperature profites from a 
const.mt temperature gradient for immersion. The pro- 
files were aligned and scaled vertically as described in 
the text. The detailed agreement is quite striking, and 
is best between CTIO and CLCO. the most distant 
stations. We attribute the poorer a^eement of the low> 
frequency components of the ESO data to known difR- 
culties in guiding during immersion. Features labeled 
A through D. discussed in the text, are proCjced by 
the correspondingly labeled features in the light curves 
shown in Fig. 9. 


that the altitude scale applies only to the 
CTIO data: the other profiles were lineariy 
stretched or compressed in altitude as dic- 
tated by the s|^e alignment. 

The AT profiles show close agreement, 
especially fw CTIO and CLCO, the two 
most distant stations. All profiles have a 
±S°K wavelike feature of similar phase and 
vertical wavelength of about two scale 
heights. The smaller-scale variations arc 
the results of spikes on the light curve, and 
there is good, but not perfect, agreement 
here as well. The ESO profile deviates most 
strongly from the others below an altitude 
of -60 km, which we attribute to the same 
partial loss of signal we regard as responsi- 
ble for the large error in the isothermal fit 
half-light time. Noise experiments confirm 
the plausibility of this suggestion. The at- 
mospheric Fresnel vertical scale is ~2 km 
in the relevant regions for these observa- 
tions, providing a resolution limit on the 
profiles in Fig. 8, since the assumptions of 
geometric optics are invalid at finer vertical 
scales. 

C, The Relationship of Light Curve Spikes 

and Temperature Profiles 

A light curve and its derived temperature 
profile contain the same information, in that 
each can be reproduced from the other, but 
the most prominent features in the light 
curve, the spikes, have only a minor effect 
on the temperature profile. The vertical 
wave in Fig. 8 is produced by slow modula- 
tion of the light curve. Since the light 
curves themselves are highly correlated, it 
is not surprising that the derived tempera- 
ture profiles are as well, but the correlation 
is not perfect. It is instructive, therefore, to 
examine in detail several light curve spikes 
and the temperature variations they pro- 
duce. We selected four such features in the 
high-resolution section of immersion data 
shown in Fig. 9. The corresponding fea- 
tures in the AT profile are labeled in Fig. 8. 

Feature A. This is an example of a fea- 
ture which is well correlated for two closely 
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Fic. 9. High-resotution segment of immersion occultation light curves. Tune is in seconds after the 
following reference times on IS August 1980 (LT): CTIO, 22^ 38" OltOI; ESO, 22^ 38" OSMO: CLCO, 
22^ 38" 0h!6S. Features labeled A through D. discussed in the text, produce the correspondingly 
labeled stnictures in the temperature profiles obtained by numerical inversion (Fig. 8). 


spaced stations (ESO and CLCO), but 
weakly correlated with the most distant sta- 
tion, CTIO. The trough between the two 
spikes corresponds to a defocusing region, 
which can be thought of as a convexity in 
the stellar wavefront as it passes through 
the Uranus atmosphere. The convexity is 
produced by a rarefied region sandwiched 
between two relatively dense regions. The 
intensity of a spike depends upon the radius 
of curvature of the effective lens, which in 
this case involves the quantity d^nldr, 
where n is the number density and r is the 
altitude. The corresponding temperature 
fluctuation is indicated on the AT profile. 
The differences among the profiles of differ- 
ent stations are very slight, but are suffi- 
cient to produce noticeable differences in 
the light curve. 

Feature B. This is a fine example of the 
gradual modification of spikes with distance 
along the limb. A large trough separates 


two pairs of spikes which are barely distinct 
in the CLCO data and are easily resolved in 
the CTIO light curve. The corresponding 
temperature maximum is quite pronounced 
in all three profiles. 

Feature C. According to the layering 
model, sharp spikes correspond to regions 
of strong focusing, which are regions in the 
atmosphere where the temperature de- 
creases sharply with height. In this case, 
two nearly identical spikes in the CLCO 
data migrate down the light curve and are 
delayed and altered when seen in the CTIO 
data. The delay of ~1!2 corresponds to an 
altitude difference of ~ 3 km over a horizon- 
tal separation of 144 km. 

Feature D. Two distinct spikes in the 
CLCO light curve coalesce in the CTIO 
data and produce the sharp ledge indicated 
in the AT profile. Once again, distinct dif- 
ferences in the light curve are barely per- 
ceptible in the temperature profile. 
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Other features. There are, c<Miversely, 
prominent features in the temperature pro- 
file which are inconspicuous in the light 
curve. The sharp feature in the ESO and 
CLCO profiles at 1 10 km below half-light in 
Fig. 8 is barely visible in the CTIO profile. 
The differences are caused by slight differ- 
ences in the depth of the “trough" in the 
light curve near time 120 sec (Fig. 2b). 

D. Horizontal Coherence and the Validity 

of the Inversion Method 

The spikes are clearly produced by co- 
herent horizontal structures whose charac- 
teristic length is well over 100 km. How- 
ever, the inversion method requires 
horizontal coherence over the effective 
path length in the atmosphere of ^2300 km 
for the results to be strictly valid. We have 
seen that the spikes correspond to small in- 
flections in the temperature profile, with a 
characteristic vertical extent of a few kilo- 
meters. and that the overall trend in the 
temperature profile is produced by low-fre- 
quency variations in the light curve. The 
computed CCFs are dominated by the 
spikes and are insensitive to these gradual, 
low-amplitude variations in intensity. Lin- 
ear extrapolation of the immersion CCF in 
Fig. 5 shows that the correlation reaches 
zero at a separation less than the required 
coherence scale for numerical inversion. 
However, this is misleading, both because 
of the uncertainty of the extrapolation and 
because this coherence refers primarily to 
the sharp features in the light curve. More 
multiple observations are required to deter- 
mine this important region of the cross-cor- 
relation, but it should be noted that agree- 
ment in the overall derived temperature 
profiles of widely spaced stations is itself 
evidence that the light curves are corre- 
lated. The unique mapping between light 
curves and temperature profiles obtained 
by stations separated by ~2000 km in the 
1977 Uranus occultation (Dunham et al., 
1980} guarantees that the low-frequency 
structure of the tight curves is similar. (This 


could be confirmed by synthesizing model 
light curves assuming an atmosi^teric re- 
fractivity profile corresponding to the low- 
est-order Fourier components of the de- 
rived temperature profiles.) In our view, 
the overall temperature profiles in Figs. 7 
and 8 accurately reflect the general charac- 
ter of the vertical structure in the Uranus 
atmosphere at this time and location on the 
planet. The snudl-scale variations in the 
profiles, perhaps distorted somewhat by 
horizontal inhomogeneities in the atmo- 
sphere, are typical of the true atmospheric 
structure. 

IV. CONCLUSIONS 

Comparison of simultaneous observa- 
tions of the 15 August 1980 occultation re- 
veals strong similarity in the upper atmo- 
spheric structure of Uranus over distances 
of more than 100 km. Cross-correlation 
analysis of the light curves themselves sug- 
gests a horizontal coherence scale of at 
least 140 km for the features producing the 
sharp spikes. This is inconsistent with 
models of isotropic turbulence, and any ac- 
ceptable anisotropic turbulence model must 
be able to account for axial ratios of 60. The 
spike correlation diminishes as the occulta- 
tion proceeds, due in part to altitude varia- 
tions of the atmospheric features and to 
slight changes in the wavefront curvature. 

The results give credence to the tempera- 
ture profiles obtained by numerical inver- 
sion, which can be used to investigate the 
structure and dynamics of the upper atmo- 
sphere of Uranus. The origin and maximum 
extent of the horizontal layers remain to be 
determined, but an upper limit to the 
strength of the quasi-isotropic component 
of turbulence can be obtained, which will 
provide a means of comparing the micro- 
structure of the upper atmosphere of 
Uranus with other planets. In a future 
publication, we will compare the power 
spectrum the intensity variations of these 
aiul other occultation observations with the 
predictions of weak-scattering turbulence 
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theory and with numerical realizations of 
strong, anisotropic turbulence. The great 
distance of Uranus from the Sun may ac- 
count in part for the relative quiescence of 
the atmosphere, as discussed by French ei 
al. (1982), and we must await further obser- 
vations before concluding that these results 
are applicable to all planetary atmospheres. 
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ABSTRACT 

Two occultation obsorvatlona of the upper atmosphere of 
Uranus are reported: (1) the 15-16 August 1980 occultation of 

KJi 12 from Cerro Tololo Inter-Amerlcan Observatory » and (11) 
the 26 April 1961 occultation of KME 13 from the Anglo- 
Austral Ian Telescope and the 1-meter telescope of the 
Australian National University. Mean atmospheric temperatures 
of 134 ± 15K for the 1980 event and 132 t 15K for the 1981 
event are derived from the light curves. A comparison of all 
available Uranus occultation data since March 1977 suggests 
that the mean atmospheric temperature of Uranus has changed 
significantly, with a typical variation of ISK yr'~^. We 
Investigate plausible energy sources that might account for 
such large temperature variations. We conclude that molecular 
and eddy diffusion, and atmospheric dynamics, are as important 
as radiation to the upper atmospheric heat balance. There is 
evidence that no significant radiative cooling had occurred at 
two sub-occultatlon points that had been in darkness for more 
than 0.5 terrestrial yr. , suggesting upper limits to the 
mixing ratios of CH4 and C2H2 of 6 x lO"^ and 5 * 10”’, 
respectively. The consistently close agreement of immersion 
and emersion mean temperatures for each occultation, in spite 
of apparently large secular changes in the atmospheric mean 
temperature, suggests that effective meridional transport 
occurs on Uranus. A continuing program of occultation 
observations in future years should reveal r^ether this 
pattern is significant. 


I. INTRODUCTION 
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Stellar occultatlon observations have revealed that the 
upper atmospheres of Jupiter (Veverka et a^« t 1974a), Uranus 
(Elliot and Dunham, 1979; Dunham ^ , 1980), and Neptune 

(Veverka et , 1974b) do not exhibit the quasi -Isothermal 
characteristic of stratospheres purely In radiative equilibrium 
but have strong temperature variations of large vertical and 
horizontal extent. By monitoring the thermal structure of 
these atmospheres, we can develop an observational context for 
a better understanding of upper atmospheric energy balance and 
dynamics. 

In this paper, we describe the results of two recent 
occultatlon observations of Uranus and investigate the 
relative importance of radiation, molecular and eddy 
diffusion, and wave dynamics to the thermal structure of the 
Uranian stratosphere. On 15-16 August 1980, Uranus occulted KM 
12, the twelfth star on the list of Klemola and Marsden (1977). 
The event was widely observed, and based on a comparison of 
light curves obtained from Cerro Tololo Interamerican 
Observatory (CTIO), European Southern Observatory (ESO) and 
Cerro Las Campanas Observatory (CLCO), we have concluded that 
the Uranian atmosphere is not strongly turbulent, and has the 
same vertical structure over large horizontal distances (French 
et al . , 1982). On 26 April 1981, Uranus occulted KME 13 
(Klemola, Mink and Elliot, 1981). We observed the event with 
the Anglo-Australian Telescope (AAT) and with the 1-meter 
telescope of the Australian National University (ANU). We 
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describe the observations in the next section, and present the 
derived atmospheric temperature profiles in $111. In $1V, we 
consider the relative importance of plausible energy sources 
to the upper atmospheric heat balance, and in the final 
section, we summarize our understanding of the thermal 
structure of the Uranian upper atmosphere. 


11. Observations and Data Reduction 


Observations of the 1S--16 August 1080 oceultation of KM 
12 by Uranus and its rings were carried out in the K band (X 
2.2 um, AX " 0.4 um) with the CTIO InSb system mounted on the 
4 m telescope (Progel and Elias. 1980). The general 
characteristics of the observations have been described in 
some detail in our analysis of the ring data (Elliot et al. , 
1981). The immersion and emersion light curves are shown in 
Figs. 1 and 2. During emersion, cirrus clouds were present 
and were responsible for the variations in signal level near 
00:16:30 UT. The immersion signal is much cleaner, and we 
shall devote most of our attention to this event. The light 
curves are characterized by numerous sharp splices. From a 
detailed comparison of observations from several stations, we 
find that the spikes are not caused by isotropic turbulence, 
but are produced by layered atmospheric structure which 
extends at least 140 km along the limb of the planet (French 
et al. . 1982). 

We observed the 26 April 1981 oceultation of KME 13 by 
Uranus and its rings with the AAT infrared system in the K 
band (X ■ 2.2 pm; Harton and Allen, 1980). The ring events 
will be decribed in a future publication. The immersion and 
emersion light curves obtained from AAT are shown in Figs. 3 
and 4; the immersion event as seen from the 1 m telescope of 
the Australian National University (ANU) is shown in Fig. 5. 
Because the ANU observations were made at 0.8 pm, the signal 
was contaminated by light from Uranus and the data are 
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significantly noisier than the AAT data, which were taken in 
the 2.2 um band, where the contribution from Uranus is snail; 

I 

however, the unocculted stellar flux was roughly 4% of the 
background at the ANU. 

On the basis of a simultaneous solution for the best 
fitting coplanar, processing elliptical orbits of the nine 
confirmed Uranian rings, using the pole of the ring plane as 
determined by Elliot et (19U1), we determined the 
occultation geometry as shown in Pig. 6. From the event 
astrometry, and assuming a planetary oblateness of 0.024 
(Elliot et al« , 1981), we computed v^, the event velocity 
perpendicular to the planetary limb. This is tabulated in 
Table I for the midpoint of the occultation, defined by ti^2* 
the time at which the best>fitting Isothermal light curve 
reached half -intensity . 


All of the light curves contain fluctuations that are 
caused partly by noise and partly by Uranus' atmosphere. To 
determine the contribution due to Uranus, we first obtained a 
measure of the rms noise in the occultation data, according 
to the relation 


c 


iMl 

n* 


1/2 


1 

N-1 


N 

I (n(tj) - 
J-1 




1/2 


( 1 ) 


where c is the rms noise in the light curve for a one-second 
integration, expressed as a fraction of n*, the counts per 
second for the unocculted star; n(tj) is the mean counting rate 
in counts per second for the integration bln of duration 
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At, N is the total number of integration bins used in the 
calculation, and fl is the mean counting rate in counts per 
second for all N integration bins. 

For the Ift-ltt August 1980 data, we computed e using the 
parameters described in French et (1982). For the AAT 

data of 20 April 1981, we used a total of 50s of data (at a 
time resolution of At = 0.1s) prior to immersion, when the 
onlv variations in the light curve would be due to noise, and 
computed the average e for a running mean of 4s duration 
(corresponding to N = 40). The results are given in Table I. 
Because emersion was observed at dawn, the post-emersion light 
curve at AAT was rising rapidly, making it difficult to 
estimate e during the occultation, and we have assigned the 
immersion AAT noise level to the emersion data as well. For 
comparison with other occultation data, we have computed a 
dimensionless measure of the noise level of an occultation, 
given by e(yx/H)^/^ (French et , 1978). This is the rms 
noise, as a fraction of the unocculted stellar flux, over the 
time interval during which the star probes one scale height in 
the atmosphere. The results are given in Table I. 

Once the contribution of noise to the total fluctuation 
strength is known, the intrinsic rms spike amplitude in the 
light curvt? due to Uranus' atmosphere can he computed, 
according to: / \ 






e^/At 


1/2 


( 2 ) 


7 


ORIGINAL PAGE IS 
OF POOR QUALITY 


where is the observed residual of the light curve from 
the local mean, normalized by the unocculted stellar intensity. 
For the 1980 observations, a data interval between 10® before 
the half light level and 20® after the half light level was 
used for immersion, and a symmetric interval for emersion. The 
corresponding interval, scaled by the corresponding v^, was 
used for the 1981 observations: between 23 S before half light 

and 45® after half light. The results are given in Table I. 
Photon noise and terrestrial scintillation are not dominant 
sources of fluctuation power, except for the ANU immersion 
data, the noisiest observation. The rms spike amplitudes are 
very similar 0.2) for all events except for 16 August 1980 
emersion, where the rms amplitude is significantly larger. 


ORIGINAL PAG£^ 

OF POOR QUALITY 

- 8 - 

III. Thermal Structure of the Uranian Upper Atmosphere 

The isothermal model light curves fit to the data in Figs. 
1-5 provide a measure of the mean temperature of the Uranian 
upper atmosphere. For the perpendicular velocities given in 
Table I and a value of 2.16 K km“^ for the ratio of temperature 
to scale height, we obtained the temperatures given in Table 
II. Since the atmosphere is manifestly non -isothermal 
is non-zero), the formal errors of the fits give no clue to the 
uncertainties of the temperature estimates. To proceed 
further, we adopted the usual assumptions of the numerical 
inversion method to obtain temperature profiles from occulta- 
tion light curves. The method is described in French ^ al . 
(1978). The critical assumption is that the atmospheric 
density varies only with radius. For Uranus, the observations 
of the 15-16 August 1980 occultation support the validity of 
this assumption, since details of the atmospheric vertical 
structure persist essentially unchanged over large horizontal 
distances (French £t ^. , 1982). We adopted an acceleration 
due to gravity of 818 cm and an atmospheric refractivity at 

STP of 1.26 X 10”*^ (corresponding to an atmosphere of 90% H 2 
and 10% He by number). The levels corresponding to full and 
zero stellar intensity were determined from the isothermal 
fits. 

The temperature profiles we obtained are shown in Fig. 7. 
The inversion procedure is particularly sensitive to noise in 
the light curve at the onset of the occultation, and because 
of the significant signal variations due to clouds during 
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emersion on 16 August 1980, we have not shown the CTIO emersion 
temperature profiles, which we regard as unreliable* To deter- 
mine the effects of noise in the data on the derived temperature 
profiles, we superimposed on an isothermal light curve actual 
noisy data prior to immersion for each occultation and obtained 
the corresponding temperature profiles. For CTIO and AAT data, 
the noise caused only low-frequency variations in the temperature 
profile of Less than 1 K amplitude. However, the ANU data 
are sufficientlv affected by terrestrial scintillation and 
background light from Uranus that the noise-induced temperature 
variations in the derived profiles are about 5 K in amplitude. 
Consequent iv , we regard these inversion results as unreliable. 

The temperature profiles in Fig. 7 are aligned vertically 
at the altitude corresponding to the time of half-light on 
each light curve; the number density at the half-light 
altitude is cm“^ and the pressure is ~ d ub. The mean 

temperatures of the atmosphere, computed over the interval 
from db kin above to loU km below the half-light altitude, are 
given in Table II. This is the region in which the inversion 
results are most reliable: at higher altitudes the initial 

condition of the inversion strongly influences the results, 
and deeper in the atmosphere, uncertainties in the baseline 
correspond ing to zero stellar intensity are important. 

because the ma.ior noise source is terrestrial scintillation 
rather than photon noise, and because the assumptions of the 
inversion method are not absolutelv satisfied, it is difficult to 
quantify the uncertainty in the overall placement of the 
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temperature profiles, but we estimate them to be accurate to 
±15K. This estimate takes into account station-to-station 
variations for multiply observed occultations, as well as 
differences resulting from Independent analysis of the same 
observations by different investigators. We emphasize that the 
variations in the temperature profiles are more reliable than 
their means or gradients and accurately represent the true 
atmospheric state to the extent that the atmosphere can be 
regarded as perfectly layered. Although the full horizontal 
extent of the stratification has not yet been determined, it 
has been demonstrated that the effects of isotropic turbulence 
in the Uranus atmosphere do not invalidate the inversion 
procedure (French £t tQ • , 1982). In what follows we assume 
that the derived profiles correspond to the true atmospheric 
structure and investigate the implications for the energy 
balance in the stratosphere. Additional multiple station 
occultation observations are necessary before more definitive 
conclusions can be drawn. 

Every temperature profile obtained to date of the Uranian 
upper atmosphere reveals temperature variations of from ±5K to 
±30K with a characteristic vertical scale of two scale heights. 
Multiple observations of the 15-18 August 1980 occultation 
reveal that these are not the artifacts of atmospheric 
turbulence, but are genuine characteristics of the layered 
atmosphere (French et rl . , 1982). The immersion and emersion 
AAT profiles have quite different structure, with the large 
amplitude immersion variations reminiscent of the 10 March 1977 
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temperature profile (cf. Pig* 5b, Dunham et al«. 1980). Also 
present are variations of a fee tenths, of a degree K over an 
altitude of a few km. These smaller-scale features are not 
due to noise, since for the 1980 data the structure is matched 
in detail in the temperature profile obtained from the CLCO 
immersion data, separated by 144 km along the limb (cf. Fig. 

7, French et aJ . , 1982). 

The upper atmospheric structure of Uranus has been 
determined from observations of four separate stellar 
occultations between March 1977 and April 1981. In Pig. 8 we 
compare the mean temperature as determined from numerical 
inversion and isothermal fits to the light curve, as reported 
for each of these events: the 10 March 1977 results are from 

Dunham e^ (1980), the 10 June 1979 results from Nicholson 
et al . (1982), and both the 15-16 August 1980 and the 26 April 
1981 results are from this work. In each case, the mean 
temperature obtained by numerical inversion was determined 
over the altitude interval from 25 km above the half-light 
level to 150 km below the half-light level. Two features are 
worth noting: first, the average temperature obtained by 
inversion varies by several tens of K over the four year 
interval for which data are available: second, the variation 
from occultation to occultation is generally larger than the 
difference between immersion and emersion mean temperatures 
for a given date. 

Are these characteristics significant, or are the 
uncertainties in the results sufficiently large that they are 
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consistent with an unchanging mean atmospheric temperature? A 
spirit of caution is in order when considering the results in 
Fig, H, It is important to bear in mind that isothermal fits 
to non-isothermal light curves may give systematically 
erroneous mean temperatures. For example, for the March 1977 
results, both immersion and emersion isothermal fit tempera- 
tures are -20 K warmer than the corresponding inversion 
temperatures. Also, it is difficult to define the mean 
temperature preciselv when the observed profiles show large 
temperature variations with height, although we have adopted a 
uniform j)rocedure for all derived temperature profiles. In 
spite of these caveats, we regard the trends in Fig. 8 as 
provocative and we snail accept them at face value for the 
purposes of discussion. As more observations become available, 
it should become clear whether these characteristics persist, 
or are fortuitous groupings of data points. 
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I V . The Upper Atmospheric Energy HudKet 

The observations In Klg. 7 anci H raise a number of 
questions: 

• vkhat is the enerRV source for the vertical temperature 
variations? 

• wfhat are the time scales for the plausible growth and 
destruction mechanisms at work? 

• Mow important are the observed small-scale temperature 
variations to the overall energy balance? 

• Are the large-scale thermal variations a dvnamical or 
a radiative phenomenon? 

• What is the cause of the apparent secular changes in 
the mean atmospheric temperature over the observed 
four year period? 

• Why are immersion and emersion mean temperatures gener- 
ally so similar for a given occultation, in spite of 
the antipodal separation of the suboccultatlon points? 

We begin our discussion of these issues by summarizing 
the characteristics of existing atmospheric models and 
com[)aring them with the observations. Next, we consider the 
energetics of radiation, diffusion, and wave motion under 
Uranian upper atmospheric conditions; then we determine the 
diurnallv and orbi tally averaged insolation on Uranus. 

Finally, we examine the observations and indicate what further 
information will be necessary to resolve the uncertainties. 

IV, A. Atmospheric Kesults 

The most detailed models of the thermal structure of the 
Uranian atmosphere have been developed bv Wallace (1 m7d, 1P80) 
and Applebv (lUMO). For the region ot interest here (number 
density cm”^; pressure -3 lih) and for plausible 
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concentrations of CH 4 , C 2 H 2 » and C 2 Hfj, the Important energy 
balance In these models involves radiation in the methane 
bands* A simple balance between heating by absorption of solar 
radiation in the V 3 (3*3 um) band and cooling by radiation in 
the V 4 (7.7 ym) band predicts a radiative equilibrium 
temperature of 125K over the region where methane is optically 
thin. This calculation assumes that LTE applies and that other 
energy sources are not significanV. The abundances of 
effective infrared radiators such as CH 4 , C 2 H 2 . and C 2 Ht> in the 
stratosphere are unknown. Uillett and Hieke (iy77) obtain an 
upper limit of < 3 k lO""®, based on their IK 

observations and an assumed temperature ofl25K. However, 
Appleby (ly«U) obtains ncH 4 * 3 x on the basis of model 

calculations, and feels that this value is not ruled out by 
Gillett and Hieke 's (iy77) observations. By varying the CH 4 , 
C 2 H 2 , and C 2 Htj abundances and taking proper account of non-LTE 
effects, Appleby has produced a variety of models with mean 
upper atmospheric temperatures ranging from 105K to 130K (cf. 
Fig. 33, Appleby (lyBO)), However, none of these models 
exhibits the relatively rapid variations of temperature with 
altitude which are observed. Clearly, some additional phvsical 


process is at work 
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IV. B. Heating and Cooling Mechantams 

In the absence of detailed Information, we begin by 
investigating In a simple way the relative Importance of 
radiation, diffusion (both eddy and molecular), and dynamics to 
upper atmospheric heating and cooling. These are the principal 
energy sources and sinks at these pressure levels. (We can 
ignore aerosol heating because the fallout rate of aerosol 
particles, governed by Epstein drag, is very large at these 
altitudes. An enormous aerosol production rate would be 
required to maintain a significant aerosol abundance.) We 
first present the functional forms of heating and cooling 
rates, and then evaluate the expression for a specific example 
of interest. 

Radiation 

The thermal energy density, AE, of a temperature 
perturbation of amplitude AT is given by 

AE “ um^CpAT (3) 

where is the mean molecular weight in gm, n is the 

atmospheric number density, and Cp is the specific heat per gm 
at constant pressure. If the vertical extent of the 
perturbation is of order the scale height. H, then Hae 
represents the thermal energy per unit surface area contained 
in the disturbance. 

The radiative heating flux, Ffp . absorbed by an optically 
thin gas within this altitude range is 


Frh - <U>fx» 


(4) 
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where <IJ> is the diurnally averaged insolation, expressed as a 
fraction of the normal incidence solar flux outside Uranus' 
atrn >spnere per unit area and wavelength interval. An upper 
limit to the equivalent width, W, of the band, adequate for 
this scaling analysis, is given by assuming that the band is 
optically thin* 

W = bnnH/L (5) 

where b is the oand strength (in ym cm“^ am“*), n is the 
fractional abundance of the absorbing gas and L is Loschmidt's 
numoer. (Under these assumptions, W is an upper limit to the 
true equivalent width, ) 

Combining fc'qs. li-b, and approximating the flux divergence 
as we can determine the characteristic time, Xj-h* heat 

the atmosphoi’e bv radiation: 

T ~ HAhi/f ) 

or 

Trh - umHCpLT/bnf \<U> 

The corresfionding radiative heating rate is 

(dT/dt)rp a bnf ;<U>/ ymHCpL (7) 

The radiative heating rate is inder)endent of the atmospheric 
density, in this approximation, and de)iends primarily on the 
abundance of the absorbing gas and the location of the 
absorption on the solar radiation curve. It is important as an 
indication of the rapidity of changes in thermal structure in a 
stratified atmosphere due to radiative imbalance. 

Analogously, a radiative cooling time for the Uranus 
atmosfJhere can be estimated. Using ttie formalism of Strohel 
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and Smith (iy73), Wasserman (1974) estimated diurnal 
temperature variations in the Jovian upper atmosphere. We 
adoyit Wasserman 's notation and apply his Kqs. (1-4) to obtain 
radiative cooling rates for CII 4 and For each constituent 


-Ahvnn(gj^/g^j)exp(-hv/kT)Fi4.oZ (8) 

wmHCp(A 

transition probabilitv, h = Planck's constant, v 
is tlie frequency associated with the transition, n is the 
mixing ratio of the radiator, and g^ are statistical 
weiglits, k is Holtznann's constant, is the probability of 

of collisional de-oxcitat ion , Z is the temperature-dependent 
collision fretpiency, and n is the total gas density. 



where A is the 


bddv Diffusion 

bulk atmospheric motion in the form of turbulence carries 
an associated energy fl»ix, Fp, which can be approximated as 

Fpj = ympngK (9) 

where g is the acceleration of gravity and K is the eddy 
diffusion coefficient (in cm^ s“M« The corresponding heating 
rate for eddy transport is as follows: 

Here, it is assumed that the laver over which F^ varies 


significantly is one scale height deep 
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Molecular Diffusion 

The molecular diffusion flux Is 

^’fn = -iiiimtiCv(dT/dz)(kT/inim^)l/2/3o (lii) 

where Cy is the specific heat at constant volume, k is 
Boltzmann's constant, T is the mean atmospheric temperature, a 
is the collision cross-section of the dominant atmospheric 
constituent, and dT/dz is the mean vertical temperature 
{gradient (Heif, iyB6), The corresponding molecular diffusion 
cooling rate is 

ydtyn, \^BonYy^iitiy / 

Here, y is the ratio of specific heats. Molecular conduction 
time scales thus decrease roughly exponential Iv with 
increasing altitude. 

Atmospheric Waves 

Typical wave dynamical energy fluxes are given bv (French 
and Gierasch, iu74): 

I'w =* Y^/( Y-1) ( 14) 

Here, w is the vertical velocity of the propagating mode, and 
the overbar denotes a time average. The dvnamical energy flux 
depends strongly on the properties of the particular 
prof>agating wave in question. The corresponding wave dvnamical 
heating rate is 

= Ywp/ ( Y-1 ) yncpH (15) 



- ly 


CM 
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An implicit assumption of Eq. 15 is that the wave deposits a 
significant fraction of its energy in the region in question. 


I V . c , Solar Energy uudget 

Armed with the preceding expressions, we can evaluate the 
relative importance of these energy sources and sinks to the 
thermal balance of the upper atmosphere. Before proceeding 
with a specific example appropriate for the occultation region, 
however, we shall provide a context tor discussion of time 
constants bv computing the diurnallv and orhitally averaged 
solar flux at Urantis as a function of latitude. 

For a planet whose rotation period is irmcn shorter than 
its orbital period, the diurnallv averaged insolation at 
latitude 0, in units of the normal-incidence solar constant at 
the planet, is given bv 

<iJ> = fzh, e, 0 p, ^)u_i(z)dX ( 16 ) 

where 

7, (X, 9, 9p, = cos i^(cos0p sinX - sin6p cosX sin 9) 

- sin4» COS0 cosX (17) 

In these expressions, 0p is the angle between the rotation axis 
and the orbital plane, X is longitude, is the angle between 
the sun and the pro.iection of the rotation axis into the 
orbital plane (pence ^ varies as the planet moves in its 
orbit), and u_i(z) is the Heaviside step function: 

( U z<U 


«i_l ( z ) 


I 


z>U 


( 1 «) 



tiinnilarlv, tho latitude dependence of the insolation averaged 
over the entire orbit can be expressed as 

<A> * /<U(0, Op, ♦)>d^ (ly) 

where, again, the average insolation is expresed as a traction 
of the normal incidence solar constant at the planet. for 
simplicity, we have assumed that the orbit is circular. 

Fig. y shows <y> for the geometry appropriate for the dates 
of observed Uranus atmosphere occul tations . The sampled 
sub-occultation latitudes are indicated bv dots. Table III 
lists <1J> for specific occultation observations. Notice that, 
because of the orientation of Uranus’ pole and the phase angle 
of the observations, two events probed regions of the atmosphere 
which had received no solar radiation for many planetarv 
rotations: the 10 June 1979 CLCO emersion sub-occultation point 

had been in darkness for O.b terrestrial yr ; the 15 August 1980 
immersion CTIU sub-occultation point for *-0.75 yr. The Uranian 
Jay is ~lt> hours long (Elliot e_t al_. , 1981). The orbitally 
averaged insolation, <A>, is shown in Fig. 10. Because of the 
extreme obllquitv of Uranus, the equatorial regions sampled by 
tho occultations receive significantly less radiation than the 
polar regions. The orbital period of Uranus is 84 vr. 


i V , I ) . Heating and Cooling Kates for Urantis* Upper Atmosphere 
We now compute energv fluxes and time constants for the 
specific circumstances appropriate for the occultation region. 


k!l 
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iVe asRumt? Applc'hv's (iy«0) mixinK ratios of ncH/j “ [L'H4]/[H2] ** 
3 X 10-5 and nc^H^ “ f 1 / ( ‘lii I = 10"®, althoiiRh these must be 
regarded as only educated guesses, and adopt the parameter 
values given in Table IV. The 0.0 ym (v^) methane hand 
strength is computed from Applebv's (1980) careful evaluation 
of available band intensities. The radiative heating rate in 
the CH4 V3 band is, from liq . (7) 



0.5 



<D> 


5 X 1U“5 


1 

0 

N 

• 

0 

1 


vr 


-1 


( 20 ) 


and the corresponding characteristic heating time is 


^rh 




T 


12oK 


0 X 10-5 
hCH4 


0.25 

< 0 > 


vr 


( 21 ) 


The cooling rates for radiation bv CH4 and C 2^2 
calculated from Eq. (8) and are shown in Figs. 11 and 12 for a 
range of atmospheric temperatures and number densities, and 
selected values of the cooling rates and corresponding radiative 
relaxation times are given in 'fable V. For temperatures and 
densities comparable to those sampled during the occultation 
(T ~ lOO-lbO K , n 10^ ^-10^5 cm“3), radiative cooling times 
are significantly longer than a vear, particularly for low 
tempf>ra til res and densiti<?s. Naturally, the greatest 
uncertainties in those figures are the minor constituent 
abundances, ncil^ ®^hd hc^H^* shall return to this point. 

The eddv iliffusion heating rate is, for the conditions of 
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the occultation and for parameter values in Table IV: 


/ dT 
dt 


2. M 


e 


K yr”^ 


lU^Ccm^s’'^ ) 

and the characteristic eddv heatinR time Is 

10 ^ ( cm^ s~ M 


( 21 ) 
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K 


vr 


( 22 ) 


For molecular diffusion, the corresponding relations are 

10^‘*cn~3 


{ 


dt/r 


m 


K vr“^ 


(23) 


and 


■ m 


10 ^ **cm” ^ 


vr 


(24) 


for a vertical temperature gradient of 0.2 K km"^. Typical 
values for these quantities are given in Table VI. At low 
nump(?r densities, the molecular diffusion time is much shorter 
than the radiative relaxation time, and even as low an eddy 
viscosity as k = 10 ^ cni^ s“^ yields a time constant shorter 
than typical radiative relaxation times at 12bK. h’vidently, 
radiative heating and cooling are slow processes in the strato- 
sphere of Uranus. The lack of detectable cooling over the long 
times during which the 1470 t.'LCO emersion and the 1980 CTIt» 
immersion sub-occul tat inn iioints were in darkness is consistent 
with the adopted values of ncjH^ ^*-2**2* however, if ncH^ 

were as large as H x 10 “^ and nCvH> were as large as b x 10“^, 
an overall stratospheric cooling of about 20K would have taken 


place over the long nighttime, and such a cooling is not 
evident in Fig. 8 . We therefore rule out such large 
abundances, unless there is effective horizontal mixing. 

Kadiative equilibrium calculations for Uranus make precise 
predictions about the thermal structure because of the strong 
temperature dependence of the 7.7 pm CH 4 emission; the Wien 
region of the blackbody curve of Uranus is extremely steep in 
this wavelength range. This acts as a highly effective 
thermostat, since radiative cooling increases rapidly as the 
mean temperature increases. However, the energy content of the 
heating radiation is low because of the great distance of 
Uranus from the sun. Consequently, even as inefficient a 
process as molecular diffusion is of comparable importance to 
the energy balance of the upper atmosphere, and we conclude 
that the atmospheric structure cannot be described adequately 
in terms of radiative equilibrium alone, regardless of the 
detailed properties of putative minor atmospheric constituents. 

Other energy sources must be responsible for the observed 
temperature variations. What are the possibilities? The 
formation and growth of a thermal perturbation is possible only 
if the energy flux divergence is non-zero. Since molecular 
diffusion always attempts to create an isothermal state, a 
lower limit to the required energy is that the energy flux 
divergence exceed that due to diffusion. The corresponding 
heating rate must be at least as large as that due to molecular 
diffusion, typically 1-lOK yr“^ in the occultation region (cf. 
Table VI). In other words, anv energy source large enough to 
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overcome the effects of molecular diffusion and produce the 
observed thermal structure must also be large enough to 
Influence significantly the basic state temperature of the 
upper atmosphere* Thus, although molecular diffusion has 
little direct effect on the basic thermal state in a radiative 
equilibrium model, the diffusion flux divergence implied by the 
observed temperature gradients gives a measure of the energy 
fluxes important to the thermal balance, and these are at least 
as important as radiative energy sources* Under these 
circumstances, it is not surprising that the observed mean 
temperature does not match the precise predictions of radiative 
equilibrium, nor that the mean atmospheric temperature is 
observed to change so quickly (Fig. 8). 

Heating by turbulent energy can be significant if the eddy 
viscosity is sufficiently large and if there is a localized 
source of energy to feed the eddies. Such a region could be 
produced by shear instabilities or by dissipation of vertically 
propagating waves, for example. The overall structure in Fig. 

7 cannot be an inertia-gravity wave, because Dunham et al * 
(19HU) have shown that such a mode would have an enormous 
energy flux. However, there are many small-scale features in 
Fig. 7. These features extend for at least 144 km along the 
limb for the lb-lb August 1980 observations. If these are 
regarded as inertia-gravity waves, then the dispersion relation 
given in French and Oierasch (1974) can be used to estimate the 
associated energy flux. For example, assuming a vertical 


wavelength, Xy, 


of b km, a horizontal wavelength of 144 km, and 
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a temperature perturbation of 0.5 K, at n the 

corresponding energy flux is ■ 0.024 erg cm"*^ 8”^, and the 
associated heating rate is (dT/dr)^ ~ 1 K yr"^. Such short 
vertical wavelength modes cannot propagate when the molecular or 
eddy diffusivity exceeds a critical value, which for this 
example Is 3 x lO** cm^ s~^. The waves then dissipate their 
energy and heat up the atmosphere. 

This is an attractive possibility because it accounts for 
both the observed fine structure of the temperature profile and 
satisfies the energy requirements. The mechanism requires a low 
eddy viscosity (K < 3 x 10** cm^ s**^) to permit waves generated 
deeper in the atmosphere to propagate into the occultation 
region without significant dissipation. (Such small values of K 
are not unheard of. For example. Tomasko (1974) determined that 
K ■ 2 X lo** cm"2 in the Jovian stratosphere at pressures 

between 20 and 500 mbar. ) Unfortunately, we have no information 
about the stratospheric eddy viscosity for Uranus, but it is 
conceivable that vertical mixing could be inhibited in the 
atmosphere of a planet which has no significnt internal heat 
source and which is so far from the sun. Stone (1973) shows 
that the troposphere of Uranus is expected to be remarkably 
stable, with very weak turbulence and a Richardson number in 
excess of 1000. Under such conditions, vertical mixing by 
eddies should be minimal. Alternatively, the temperature 
structure may reflect dissipation of turbulent energy resident 
in a stratospheric region of strong shear, but little further 


can be said in the absence of more information 
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One of the intriguing aspects of Pig, 8 is the generally 
close agreement of immersion an^l emersion mean temperatures, 
compared to the larger variation in temperature with time. If 
these observations are representative, then they suggest that 
the atmospheric circulation is sufficiently strong to 
isothermalize the stratosphere, in spite of the large 
differences in dlurnally averaged insolation (Pig. 9). 

Although radiative time constants are long in the occultation 
region, they are shorter deeper in the atmosphere, and a 
typical radiative relaxation time at the 2 mb level is seven 
days (Appleby, 1981, private communication). Consequently, the 
variations in <D> shown in Pig. 9 may be important drives 
to circulation at deeper levels which ultimately influences 
the upper atmospheric thermal structure on time scales shorter 
than the local radiative time constant. It is unclear what the 
global mean circulation is on a planet with a deep atmosphere 
whose poles receive more energy than the equator (Fig. 10), and 
even less clear how important these tropospheric dynamics are 
to the stratospheric structure. 8tone (1976) suggests that the 
deeper layers where solar radiation is absorbed will have a 
moderating effect on the higher layers where the time constant 
is considerably shorter. 

Many of these issues will remain unsettled until Voyager 
encounters Uranus in 1986, but some of them may be clarified 
considerably in earlv 1982, when three high quality atmospheric 
occultations of Uranus will he observed within six weeks 
(Klemola, Mink and Klliot, 1981). Prom these observations, we 


can learn whether the agreement between immersion and emersion 
temperatures persists, and can chart more accurately the 
secular changes in the mean atmospheric temperature. Multiple 
observations of each event will permit the horizontal scales of 
the temperature perturbations to be determined, allowing more 
precise identification of the possible dynamical modes they 


represent 
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V. CONCLUSIONS 

Observations of both the mean stratospheric temperature 
and the detailed thermal structure of Uranus* atmosphere reveal 
that non-radlatlve sources of enerKV Are important to the heat 
balance, with heating rates of tens of K yr**^. Kadiative time 
constants are longer than the time scales over which signifi* 
cant changes in the mean stratospheric temperature have been 
observed* The possibility that the upper atmosphere is 
Isothermallzed over time scales shorter than a few months 
suggests that the general circulation may be very effective at 
meridional transport of energy at these levels. Under these 
circumstances, discussions of a constant amplitude temperature 
Inversion of the upper atmosphere are simplistic, and ignore 
important, and perhaps dominant, non-radiat Ive energy sources. 
The complexities of the picture can be sorted out only when 
more precise measurements of the abundances of radiatively 
Important minor constituents such as C2H2 and C2^ti have been 
made, and when additional occultatlon observations confirm or 
deny the existence of the patterns presented by the 
observations available so far. At that time, it will be 
appropriate to construct more realistic atmospheric models 
which take proper account of these non-radiative effects, in 
order to determine the sensitivity of the thermal structure to 
changes in thermal and dynamical forcing. 
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TABLE X 


Ll^t CttTV Propertiao and Occultation Timing 


Date 

Station 

Event* 

vj^(km-s“^)^ 

C<s^/^) 

e( vj^/H) 


ti/2(OT) 

15 Aug 1980 

CTIO 

1 

6.811 

0.007 

0.002 

0. 18 

22t38t09.07 

16 Aug 1980 

CTIO 

E 

6.888 

0.008 

0.003 

0.28 

00t16i26.42 

26 April 1981 

ANU 

I 

15.248 

0.123 

0.058 

0.18 

19t33t07.80 


AAT 

I 

15.248 

0.007 

0.003 

0.17 

19:33t07t12 


AAT 

E 

15.196 

(0.007)‘* 

(0.003) 

0.14 

20:04(43.85 


Notes 

* X ■ Immersion, E ■ Emersion 

^ Event velocities confuted for an assumed planetary elliptlcity of 0.024. 

^ Scale height of 70 kra assumed. 

e for AAT emersion was assumed to be the same as for AAT imnersion, since 
sunrise prevented an independent determination of e from post**emersion 
data. 
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TAbLI:; II 

Mean Atmospheric Temperatures 


Date 


Station 

Event*^ 

T(k) 

(Isothermal Fit) 

T(K)b 

( Inversion ) 

15 August 

1080 

CTIO 

I 

151 

156 

10 August 

1980 

CTIO 

E 

156 

c 

26 April 

1981 

ANU 

1 

(121)d 

(119)<1 



AAT 

1 

129 

138 



AAT 

E 

131 

129 


Notes 

^ 1 = immersion, E » emersion. 

b The mean temperatures obtained from numerical inversion 
were averaged over an altitude interval of ~175 km, from 
175 km, from 25 km above to 150 km below the half-light 
level. The midpoint of this range is near number 
density 3 x cm"^. 

c Inversion of this emersion event was not reliable 
because of clouds. 

^ These data are significantly noisier than the AAT 
observations, (cf. Table I). 



TABLE III 


Diurnally and Annually Averaged Insolation at the 
Obaerved Sub-Occultatlon Latltudaa on Uranus 


Date 

Station 

Event 

Sub-Occultation 
Latitude on Uranus 

<D>^ 

<A>* 

10 March 1977 

KAO 

I 

4t7 

0.23 

0.21 


KAO 

E 

-18. 1 

0.08 

0.22 


CAPE 

I 

7.4 

0.25 

0.21 


CAPE 

E 

-20.8 

0.06 

0.23 

10 June 1979 

CLCO 

I 

26.2 

0.39 

0.24 


CLCO 

E 

-30.4 

0.0* 

0.25 

15-16 August 1980 

CTIO 

I 

-26.3 

0.0* 

0.24 


CTIO 

E 

17.9 

0.30 

0.22 

26 April 1981 

AAT 

1 

8.4 

0.19 

0.21 


AAT 

E 

-18.6 

0.004 

0.23 

Notes 






^ <D> is the 

diurnally 

averaged 

insolation in units of 

the normal 



incidence solar constant at Uranus, for the observed sub-occultation 
latitude. 

* The observed sub-occultation latitude was in darkness throughout a 
full rotation of the planet. See text for details. 

* <A> is the orbitally averaged insolation in units of the normal 
incidence solar constant at Uranus, for the observed sub-occultation 
latitude. 
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TABLK IV 


Adopted Parameter Values 


umn 

Cp 

Y 

b 

(3.3 um CH 4 band) 

nCijHij 

B 

o(H2 ) 

dT/dz 

fx 

(at X * 3.3 ym) 

T 


3.7 X 10-2** g 
1.17 X 10® erg g“l K 
1.47 


0.25 pm cm”^ am“^ 
3 X 10“® 

2 X 10-® 

BIH cm s“^ 

6.9 X 10“^® cm2 
0.2 K km“^ 

54.2 erg-cm“2 g-l 
125 K 



Radiatlva Cooling Ratea and Relaxation Timca 
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TABLK VI 

Thermal Relaxation Times* 


Process 

n(cm“^ ) 

K(cm2 s“l) 

d'f/dt(K yr“l) 

T(yr) 

Molecular 

10^2 

.. 

3t>0 

0.45 

Diffusion 

ID 13 

— 

36 

4.5 


lOl- 


3.6 

45 


lol 5 


0.36 

450 


Kddy 

— 

lO** 

.38 

330 

Diffusion 

— 

10 5 

3.8 

3.3 


— 

10^ 

380 

0.33 



107 

3.8 X 103 

3.3 X 10-2 


— — 

108 

3.8 X lO** 

3.3 X 10-3 


* A temperature of 125K is assiimed 



KIGUKK CAPTIONS 


Pig. 1. Immersion light curve of the 15-lH August 19S0 

occultation observed from CTIO. The light curve begins 
at 'A2 : [id : 4i) » 0 (UT) and is plotted nt 0?ii resolution. 

The best fitting isothermal light curve is also shown, 
corresponding to a temperature of 151K. 

Pig. 2, Emersion light curve of the lb-lb August 19»0 

occultation observed from CTIO. The light curve begins 
at 00:12:30.0 (UT) 16 August lOHO and is plotted at 0?2 
resolution. The best fitting isothermal light curve is 
also shown, corresponding to a temperature of 156K, The 
strong undulations at the end of the plotted data were 
caused bv clouds. 

Pig. 3. Immersion light curve of the 26 April 19H1 

occultation observed from AAT. The light ctirve begins at 
19:31:44.4 (UT) and is plotted at o53 resolution. The 
best fitting isothermal light curve is also shown, 
corresponding to a temperature of 129K. 

Pig. 4. Lnersion light curve of the 26 April 1981 occultation 
observed from the AAT. The light curve begins at 
2 (J:uu: 29.4 (UT) and is plotted at «)?3 resolution. The 
best fitting isothermal light curve is also shown, 
corresponding to a t<*mperature of 131K. 

Pig. b. Immersion light curve of the 26 April 19H1 

occultation observed from ANU. The light curve begins at 


. *.*k 
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19:31:44 (UT) and is plotted at 1® resolution. The best 
fitting isothermal light curve is also shown, 
corresponding to a temperature of 121K. The data are 
considerably noisier than those shown in Fig. 4 because 
a smaller telescope was used, and the observations were 
not carried out in the 2.2 ym methane band. 

Fig. b. Occultation Geometry. The tracks of the 15 August 

1980 and 26 April 1981 events are shown for CTIO and AAT, 
respectively. The labelled dots indicate the location of 
the Uranus north pole for the two occultations. The 
immersion and emersion 1980 events probed latitudes -26* 
and ■►18“, respectivelv , on Uranus. The immersion and 
emersion latitudes for the 1981 events are ■►8* and -19*, 
respectively . 

Fig. 7. Temperature profiles of the Uranian stratosphere 

obtained bv numerical inversion of the light curves in 
Figs. 2, 4 and 5. The profiles are aligned in altitude 
at the point corresponding to the half-light time 
obtained from an Isothermal fit to the data. This 
corresponds roughly to a number density of 3 x 10^** cm"^ 
and a pressure of 6 yb. Because of the uncertainties in 
the initial condition and the choice of baselines used 
for the Inversion, the overall placement of the 
temperature profiles are uncertain by about 15K. 

Fig. 8. Occultation determinations of mean upper atmospheric 
temperatures for Uranus. The circles mark the 
temperatures obtained from isothermal fits to the 
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occultation light curves obtained on the following dates: 
10 March IH77 from the Kulper Airborne Observatory, 10 
June 1079 from Cerro Las Campanas Observatory, 15-16 
August from CTIO, and 26 April 1981 from the AAT. The 
squares mark the mean temperatures obtained from 
numerical inversion of the same data. Open symbols 
correspond to immersion and closed symbols correspond to 
emersion events. The variations in mean temperature with 
time are generally larger than the variations in 
temperature with location in the planet (immersion and 
emersion) for a given occultation. Typical uncertainties 
in mean temperatures are of order il5K. This estimate 
takes into account differences resulting from independent 
analysis of the same observations by different 
investigators, as well as station-to-station variations 
for multiply observed events. 

ig. 0. Diurnally averaged insolation on Uranus, as a 
function of latitude, for the dates of the observed 
atmosphere occultation. The ordinate is in units of the 
normal incidence solar constant at Uranus. The large 
dots correspond to the Uranus latitudes probed for each 
date, 

ig. 10. Orbitallv averaged insolation on Uranus for one 
Uranian year, as a function of latitude. A circular 
orbit was assumed. The ordinate is in units of the 
normal Incidence solar constant at Uranus. The polar 
regions receive, on average, more solar energy per orbit 
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than the equatorial regions. 

Fig. 11. Cooling rates for emmission In the 7.7 ym CH4 band, 
as a function of total atmospheric number density and 
temperature. The assumed CH4 mixing ratio is TICH4 ■ 

3 X 10-^. 

Fig. 13. Cooling rates for emission in the 13.7 um C2H2 band, 
as a function of total atmospheric number density and 
temperature. The assumed C2H2 mixing ratio is nc2H2 * 

2 X lU"®. 
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I . Introduction 

Stellar occultations are sensitive probes ot distant 
planetary atmospheres, anci n some cases, provide the only 
information we have about the structure of these remote upper 
atmospheres. Unfortunately, it has not proven easy to extract 
definitive results from the observations. The spiky structure 
of the light curves has been variously interpreted as 
scintillation due to isotropic turbulence (Jokipii and Hubbard, 
1977; Texas-Arizona Occultation Group, 1977; Hubbard, 1979a) 
and as evidence for layered atmospheric structure (Kovalevsky 
and Link, 1969; French and Elliot, 1979; French et ^. , 1982b). 
A large (and occasionally acrimonious) body of literature has 
been generated by champions of each view, but there has been 
little common ground for discussion. Proponents of turbulent 
scintillation argue that the assumption that the atmosphere is 
perfectly layered is naive and unrealistic, while their critics 
contend that the assumptions of weak scintillation theory are 
violated by the very observations the theory is used to 
explain. 

A new approach is called for: the investigation of the 

properties of light curves for atmospheres with more realistic 
properties than those considered so far. In SI I we relax the 
restrictions of weak scintillation theory and construct model 
light curves for atmospheres with strong, ar,isotropic 
turbulence using wave optics. We compare the results with 
available weak scintillation theory and with recent occultation 
data. Then, in SHI, we examine the effects oi large-scale 


atmospheric waves with realistic horizontal structure. We 
assess the reliability of the numerical inversion method in 
retrieving the true atmospheric vertical structure under 
circumstances of strong ray crossing and horizontal 
inhomogeneities. Next, we compare the nature of model light 
curve spikes generated using geometric optics and wave optics. 
In the last section, we present our conclusions and indicate 
fruitful areas for further work. In an appendix, we provide a 
list of definitions of the symbols employed in this paper. 

1 1 . Strong Scintillation in Stellar Qccultations 

The most striking feature of typical stellar occultation 
light curves is the presence of numerous sharp and intense 
spikes, produced by non-isothermal structure in the planetary 
atmosphere. Young (1976) first articulated the view that 
spikes were the manifestation of turbulence, and subsequently 
the theory of turbulent occultations has been well developed by 
Hubbard, Jokipii, and Wilkening (1978, henceforth referred to 
as HJW), and by Haugstad (1978a, b,c and 1979). and Eshleman 
and Haugstad (1979). [Although there has been a long-standing 
debate between these two groups concerning certain details of 
the theory, (Hubbard, 1979; Eshleman and Haugstad, 1979, and 
references cited therein) these differences are not important 
for the present discussion.] Woo ^ al. (1975, 1980b) have 
developed an analogous theory for radio occultations. 

These authors restrict their treatment to weak 
scintillation, because the more general case is analytically 
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Intractable. This is a serious limitation, since the observed 
rms spike amplitude is not small compared to the mean intensity 
in most occultations , contrary to the fundamental assumption of 
weak scintillation theory. In this section, we investigate the 
properties o> strong scintillation in stellar occultations by 
generating model light curves for strongly turbulent 
atmospheres. We confirm the results of weak scintillation 
theory, and summarize the qualitative differences apparent as 
the turbulence strength in increased. Finally, we discuss the 
applicability of these results to actual occultation 
observations. 

A. Thin Screen Model 

The general problem of stellar occultation theory is to 
relate the refractive properties of a planetary atmosphere to 
the characteristic features of the corresponding light curve. 
The occultation geometry is shown in Fig. 1, where an incident 
plane wave from a star is distorted by the planetary atmosphere 
on its path toward the earth. The effect of the atmosphere can 
be represented by a thin phase-changing screen which changes 
the phase of the incident plane wave but not its amplitude. 

The thin screen approximation is well satisfied for stellar 
occultations, as HJW demonstrate in some detail. 

We adopt the coordinate system showwon in Fig 1. The star 
is considered to be monochromatic with vacuum wavelength X. 

Let ♦(x,y) represe *■. the phase delay suffered by the emergent 
ray passing through the point x,y in a reference plane 
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perpendicular to the line of sight. In the absence of the 
planet's atmosphere we would have 4 « 0. Thus: 

♦(x.y) = / dz v(x,y) (1) 

where v Is the atmospheric ref ractlvlty . The lateral 
deflection of the ray Is negligible within the atmosphere. 

The atmosphere Is assumed to be such that the refractlvlty 
can be decomposed as 

= VA(r) + vt(x,y,z) (2) 

where Is the refractlvlty of the average atmosphere (with r 
the distance from the center of the planet), and where Is 
the variation of the refractlvlty from its average value. The 
variations may represent turbulence, density waves, or other 
phenomena in the atmosphere. For the consideration of 
turbulence, we adopt an isothermal average atmosphere: 

vaCt) = Voexp[-(r-ro)/H] (3) 

where Tq is a reference altitude, Vq is the refractlvlty at 
this level, and H is the atmospheric scale height. Letting r <5 
be the altitude corresponding to the half-intensity point of 
the average light curve, we have the half-light refractlvlty 

h3/2 

(2TfR)l/2D 

where R is the planetary radius and D is the distance from 
the earth to the planet. 

Furthermore, we assume that the rms turbulent refractlvlty 
is proportional to the local mean refractlvlty throughout the 


atmosphere : 
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(5) 

It follows from the decomposition of Eq. 2 that one has from 
Eq. 1; 

♦ » ( 6 ) 
where ♦a is the phase shift due to the average atmosphere, and 
that due to the variations. For the "half-intensity” ray, 

~ 2ir(H/rf)2, where the Fresnel scale is given by rf * 
(XD)^/2. Generally, the point source wave amplitude at the 
reference screen is simply 

Es(x,y) = exp(i*A i*t) (7) 

The Fresnel small-angle approximation to the wave propagator is 
applicable, and, therefore, the point source wave amplitude at 
the earth is given by 

E(x,y) = jr/dx'dy'Es(x',y') exp|^ [(x-x ' ) 2 + (y-y')2]|(8) 

The wave intensity at the earth is 

I(x,y) = lE(x,y)l2 (9) 

For a finite angular size source or star, in which different 
parts of its surface radiate incoherently, the intensity at the 
earth is given by 

lfs(x»y) “ / Jd ®y^ ( ®x » ®y ) ^ ®x^ » y~0yD) (10) 

where S(9x,0y) is the brightness distribution of the star and 9 
is a measure of subtended angle as viewed from the earth. We 
have in mind the common situation where the star is small in 
the sense that d* = 0*D<<H, where 0>*> is the star's angular 
diameter. 

Our approach is to evaluate Eqs. 8-10 numerically to 
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produce simulated light waves for cases of Interest, and to 
determine the statistical properties of the fluctuations. The 
basic properties of the resultant curves are governed by the 

A A 

non-dlmenslonal scale factors H » H/rf ar'd d^^ » do/rf, and by 
the nature and strength of the turbulence field. Represent- 
ative scale factors for recent stellar occultatlons are given 
in Table I. 

The central idea behind the simulations is that Eqs. 8 
and 10 are of the form of convolutions, so that the fast 
Fourier transforms (FFTs) of these equations become simple 
products of the FFTs of the separate factors. In practice, the 
light curve is computed in overlapping pieces which are 
appropriately scaled and fitted together. Fig. 1 shows 
schematically how two adjacent portions of the light curve, A 
and B, are computed from overlapping segments of the upstream 
phase screen. Care is taken to insure that there is no 
aliasing during the computation, and the coordinates are scaled 
by the effective Fresnel scale in such a way that there is no 
loss of precision deep in the occultation where becomes very 
large. 

Once the point source light curve has been obtained, the 
influence of the finite source is taken into account by 
convolving the light curve with the intensity distribution of a 
uniform stellar disk. Again, this is accomplished with FFT's 
and application of the convolution theorem. 
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B. Anisotropic Kolmogorov Turbulence 

An aspect of Eqs. 7-10 which is poorly known is the 
characterization of the refractivity variations v^, or the 
phase variations ♦t* reference plane. The fact that 

observed light curves of finite size stars have prominent 
spikes provides a broad constraint on the nature of 
Furthermore, there is evidence, in some cases, of substantial 
anisotropy of whatever atmospheric structures are producing the 
spikes. As early as I960, Kovalevsky and Link proposed that 
the spikes in Neptune occultation light curves were caused by 
layered structure rather than by turbulent scintillation. More 
recent Uranus observations have provided unequivocal evidence 
of homogeneous fine scale horizontal structure extending for 
hundreds of km along the planetary limb (French ^ , 1982b). 

Woo ^ al (1980a) concluded that the turbules in the atmosphere 
of Venus have characteristic aspect ratios of 20 or more. 
Terrestrial observations are consistent with Kolmogorov 
turbulence for horizontal scales of up to a dozen km (Armand et 
, 1975). 

These observations provide the motivation for investigating 
in some detail :;he nature of strong scintillations produced by 
anisotropic turbulence-. As an extreme case of such models, we 
will consider perfectly layered Kolmogorov turbulence in which 
the refractivity varies only with radial distance from the 
planet's center, r. (We are interested for the moment in the 
characteristics of the resulting light curves, and will defer a 
discussion of the actual realization of such turbulence fields 
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ORIG’TxAL ?A^.G IS 

OF POOR QUALITY 

In real atmospheres.) Letting Pv(4x) wavenumber power 

spectrum of the refractlvity fluctuations, with as the 
vertical wavenumber, one can readily show that wavenumber 
spectrum of the phase on the reference screen Is given by 

Ps(Qx) “ C2ir)3-^ (H) 

^ I Ax i 

(Questions for R.L. : Is this valid for the geometry in 

question ?) 

For one-dimensional Kolmogorov turbulence, the refractlvity 
fluctuations have the. power spectrum 

Pv(Ax) - (Ax^ + 1 /Lc2)“XX (12) 

where Lg is the outer scale of the turbulence, assumed to be 
equal to H. 

Representative results are shown in Fig, 2 for point 
source light curves, (d = 0) with one-dimensional Kolmogorov 

A 

turbulence and H » 30, a typical value for recent occultations 
(cf. Table I). The turbulence strength is defined by A, the 
rms phase shift in radians at the altitude corresponding to the 
half-light ray. Note that from Eqs. 3-5, 

// y ^ ^ = _A ^ 

NV^'a; / 2ir H 

and 

- A rf^ exp[-(r-rn/H] (14) 

(2ir) 3/2(RH) 

The curves in Fig. 3 range from A * 0.0, equivalent to an 
isothermal atmosphere, to A = 4.0, corresponding to very strong 
turbulence. In each case, the same random phase front was 
used, scaled by the turbulence amplitude, A, so that the 
modification of individual light curve features can be observed 
as A increases. The curves are characterized by many narrow, 


(RL-Hfh 

should 

the 

expone 
be for 
this 
case?) 
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Intense spikes. 

The effects of convolution with a finite source are seen 
in Fig. 3, for cases of both weak (A <■ 0.5) and strong 
(A > 4.0) scintillation, and for stellar diameters of d* ■ 0, 
1, 2, 5 and 10 rf. Even for very strong scintillation, the 
spikes are strongly suppressed by the stellar smearing. 

(Further remarks?) 

Statistics of Spike Intensity 

A statistical measure of the scintillation strength is 
given by <(A^)^>^/2, the rms deviation of the light curve from 
the mean, isothermal light curve. Weak scintillation 
corresponds to the case where the fractional Intensity 
fluctuations at a given point in the light curve are small 
compared with unity. Equivalently, the modulation index, m, is 
small compared to unity. Here: 

" - {{My) 

where is the intensity of the mean (isothermal) light curve. 
Fig. 4 displays the rms amplitude and the modulation index of 
the fluctuations in the light curves of Fig. 2. In each box, 
the triangles correspond to the computed values, averaged over 
a distance of 2.5 H along the light curve. (Adjacent points 
are separated by 0.5 H. ) The solid curves correspond to the 
predictions of HJW [their Eqs. 36-39] for layered Kolmogorov 
turbulence, based on weak scintillation theory, and properly 
scaled for the appropriate turbulence amplitude, A. Notice 



carefully that the vertical axis scales differ from plot to 
plot. Also shown in each case is an isothermal light curve. 

Several points are worth noting. There is excellent 
agreement between the HJW curves and the simulations for A » 
0.25 and 0.50, but the agreement becomes poorer as the 
turbulence amplitude is increased. (The apparent systematic 
offset of the rms amplitude model curve and simulation results 
between ^ = 0.9 and ^ = 0.7 is due to the large (2.5 H) bin 
size employed to compute the average fluctuations, which biases 
the results where the light curve is steep. ) Similar results 
were found using turbulence spectra generated from independent 
sets of random noise, providing assurance that the results in 
Fig. 4 are truly representative of anisotropic Kolmogorov 
turbulence. 

As the turbulence amplitude increases, the modulation 
index saturates at unity. This result "has been more or less 
accepted as an empirical fact in the field of radio waves" 
(Rumsey, 1975), and has been demonstrated numerically for 
strong, isotropic turbulence in a homogeneous background 
(Marians, 1975). According to Fante (1980), "it is now 
generally accepted that in very strong turbulence, the 
normalized variance of the irradiance fluctuations approaches 
unity." This is the first demonstration of such saturation 
under stellar occultation conditions, with a strongly 
defocussing background atmosphere. 
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Single Frequency Autocorrelation Functions 

Under conditions of weak scintillation, the width of 
typical spikes increases along the light curve. This is 
because the Fresnel radius scales as in the atmosphere 

(HJW) and adjacent rays diverge as with the result that 

the characteristic spike width increases as The 

autocorrelation function (ACF) of a light curve provides a 
measure of typical spike widths, and we co? puted ACF's of point 
source light curves fpr both weak (A ■ 0.5) and strong (A 
4.0) scintillation according to the relation: 

Cd.j.T) . /^AH(f )M.j(x-t')dt' (IB) 

Which represents the cross-correlation between residuals, A^, 
of light curves i and j , at a lag, t, computed over a centered 
interval of length 2t. For the ACF calculations, i equals j, 
and t corresponds to 2.6 H along the light curve. The results 
are shown in Figs. 5 and 6. Note that the width of the ACF 
increases as as predicted, for the weak scintillation 

case, but the FWHM of the ACF reaches a maximum and then 
decreases with decreasing stellar intensity. (Further 
comments?) [Discuss ACF's of finite sources?) 

Two Frequency Point Source Correlation Fanctlons 

A measure of the importance of diffraction in the 
formation of an occultation light curve is provided by the 
cross-correlation function (CCF) of observations observed 
simultaneously at two different frequencies. In a non- 


dispersive layered medium, goemetric optics dictates that the 
two light curves would be identical, and the peak CCF would be 
unity. On the other hand, the Fresnel radius varies with 
frequency, so that the fine-scale refractivity variations are 
filtered differently in the two cases. As a result, the 
maximum of the CCF is expected to be less than unity, on the 
basis of wave optics. 

We have determined the peak CCi’ for two point source light 
curves constructed from the same turbulence field but differing 
by a factor of two in wavelength, as a function of turbulence 
amplitude and distance along the light curve. The results are 
shown in Fig. 7. The turbulence amplitude shown corresponds to 
the shorter wavelength light curve; the corresponding values 
for the longer wavelength light curves are a factor of two 
lower. HJW predict a maximum CCF of 0.63 for this case, 
independent of using weak scintillation theory. Our results 
indicate that the maximum cross-correlation does have a weak 
dependence on and that the maximum decreases as the 
turbulence strength is increased. [Further comments] 

These results apply only to the special case of a point 
source. In fact, for extended sources, the CCF maximum will be 
larger than shown in Fig. 7, as long as FWHM of the point 
source light curve is smaller than the FWHM of the stellar 
image brightness distribution. For weak scintillations, the 
effect of image convolution will eventually become unimportant, 
since the characteristic spike width increases as and 

will, at some value of 4, greatly exceed the image size. As a 
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rough guide, this condition is satisfied when ^ < d«~^. When 
scintillations are strong, the situation is more complicated: 
the spike width does not increase indefinitely as ^ decreases 
(Figs. 5 and 6), while at the same time, the two-frequency CCF 
peak is depressed as ^ decreases. In this case, the image 
convolution has an important effect deep into the occultation, 
but the Intrinsic correlation between the two point source 
light curves is low. 

The significance of these results is that weak 
scintillation theory provides only a rough limit to the degree 
of similarity of two-frequency observations. If stellar 
occultations are evidence of strong turbulence, then the 
maximum cross-correlation of the two-frequency light curves 
should be less than that predicted by HJW, when the effects of 
a finite scarce are not dominant. In practical cases, it will 
be important to take care to include such convolution effects, 
which could be simulated by generating model light curves with 
parameters appropriate for the specific occultation in 
question. 

Power Spectra of Light Curve Residuals 

Although a power spectrum is just the Fourier transform of 
the ACF, and hence contains no additional information, it does 
display more clearly than the ACF the relative importance of 
high-frequency components of light curve residuals. In radio 
occultations, power spectra have been used to determine the 
scintillation characteristics in Venus (Woo et al., 19b0a) 
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but they have not been used In stellar occultations because of 
the importance of noise and the consequent limited dynamic 
range of the usable region of the power spectrum. Nevertheless, 
it is instructive to translate the ACF results in Pigs. 5 and 6 
into the parlance of power spectra, and to investigate the 
effects of finite image size in the presence of strong 
scintillation. 

To begin, we present a theoretical power spectrum for weak 
anisotropic turbulence in Fig. tt, computed by evaluating the 
Fourier transform of HJW's Eq. 67 for the one-frequency point 
source ACF. For this case, we used an anisotropy parameter of 
p * 50 (the axial ratio of the turbules), 4 - 0.1, and H * 30; 
wavenumber one corresponds to a distance of 5.12 H along the 
light curve. The maximum of the power spectrum occurs near 
wavenumber 25, corresponding to a wavelength of about 6 rj. 

Note from Fig. 5 that the weak scintillation ACF for 4 ■ 1 is 
approximately a cosine wave of this wavelength, as one would 

A 

expect. The fringes in Pig. 8 are closely spaced because H 
<< 1, and the fringes are deep because p >> 1. 

In passing, we note that we also computed power spectra 
for weak isotropic turbulence, and confirmed that the results 
of riJW and of Haugstad (1978) are identical in this case, in 
spite of their dissimilar outward form [Eq. 67 in HJW and Eq. 
3.3 in Haugstad, 1979.) Haugstad does not account for 
anisotropy in his presentation, but he includes several figures 

A 

of intensity power spectra for H - 3 as a function of 4 . 

For comparison, we computed the power spectra of ligtt 
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curve residuals for weak (A « 0.5) and strong (A « 4.0) 

A 

scintillations, wltb finite sources of d ■ 0.0, 2.0, and 10.0. 
The results are shown in Pigs. 9 and 10. Again, waveuumber one 
corresponds to 5.12 H along the light curve. The theoretical 
power spectrum (Fig. 8) represents the average of a large 
statistical ensemble, while Pigs. 9 and 10 correspond to single 
samples from such an ensemble, and are therefore somewhat 
noisy. However, it should be noted that each occultation 
observation is itself such a sample, and that Interpretation of 
the power spectrum of actual data is subject to the same 
limitations. The weak scintillation point source power 
spectrum (Fig. 9) peaks near wavenumber 25 and has the same 
slope at higher wavenumbers as the theoretical spectrum (Fig. 
8), but the deep fringes have been averaged over and are lost 
in the noise. As the stellar image size is increased, the high 
frequency tail of the spectrum is truncated and the fringes are 
reduced, but there remains a significant peak in power near 
wavenumber 25, where the point source spectrum is at maximum. 

In contrast, the strong scintillation point source power 
spectrum (Fig. 10) does not turn down sharply until wavenumber 
100, corresponding to a wavelength of 1.5rf. Again, this is 
consistent with the results of Figs. 5 and 6, which show that 
t!ie ACF is narrower for strong than for weak scintillation, and 
consequently has relatively strong high frequency Fourier 
components. [Further comments? - Effects of changing spectral 
index? ] 
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Comparisoo Wltft Observations 

To what extent do occultatlon observations resemble the 
anisotropic turbulence simulations computed above? For 
purposes of discussion, we will consider observations of the 
1980 August IS occultatlon of KM 12 from Cerro Tololo 
Interamerican Observatory. A detailed discussion of these 
observations is given by French et (1982a), and from a 
comparison of light curves obtained by widely separated 
observers (French et al . , 1982b) concluded that the data are 
fundamentally inconsistent with isotropic turbulence theory. 

It remains to consider whether these layers, which are coherent 
over distances in excess of 60 rf, are turbulent. 

As Crane (1980) points out, the scaling arguments from 
which the Kolmogorov spectrum is derived apply to the inertial- 
convective subrange, and, for these scales, the turbulent 
fluctuations are expected to be isotropic. At larger length 
scales, in the buoyancy subrange, the turbulence is 
anisotropic, but the spectrum need not be Kolmogorov, and 
depends both on the mean atmospheric flow and the strength of 
the short-period gravity waves present in the flow (Crane, 
1980). In the absence of such detailed information about the 
upper atmosphere of Uranus, we will content ourselves for the 
moment with a comparison of the statistical properties ox the 
data and the anisotropic turbulence simulations given above, 
with an outer scale of one scale height. 

In Fig. 11, the CTIG immersior. light curve of the 1980 
event is shown, together with the rms amplitude and modulation 
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index of the fluctuations about the mean curve, computed as 

A 

before for Fig. 4. For this event, d* ■ 1 (French et al. . 
19U2a), and simulations of weak and strong scintillations for 

A 

d* • 1 are shown in Fig. 3. Of the simulations shown, the 

A 

observations most closely resemble the A ■ 4.0, d* ■ 1.0, 
although the spikes are sparser in the data than in the model 
curves. Statistical analysis of the spikes shows that the 
modulation index of the observations saturates near unity, and 
a broad maximum in the rms spike amplitude Just below 
half-light is similar to the strong scintillation results 
(A « 4.0) in Fig. 4 in shape and amplitude. It should be noted 
that the modulation index is insensitive to the low-frequency 
variations of the mean light curve that are reflected as 
large-scale vertical temperature variations in the atmospheric 
profile obtained by inversion. 

The characteristic FWHM of the ACF of the observations is 
roughly 4.1 km, or 1.7 rf (French ^ a^« , 19«2b). In contrast, 
the maximum ACF FWHM for model strong turbulence light curves 

A 

with d* » 1 is a 0.9 rf. Thus, the spikes are significantly 
broader than predicted by turbulence theory. There is no 
evidence for spike broadening deep In the occultation, in 
contrast to the predictions of weak scintillation theory. 

In summary, the observed modulation index resembles the 
turbulence simulations, but the ACF is significantly broader 
than the diffraction models predict, suggesting that in some 
respects, geometric optics may be adequate to account for the 
observations. Observations at two (or nwre) widely separated 
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frequencies are necessary before more definitive conclusions 
can be drawn. In the meantime » it is worth estimating in a 
simple way the strength of refractlvity variations required to 
produce strong scintillations such as shown in Fig. 2. 

Scintillation Strength and Atmospheric Irregularities 

For isotropic Kolmogorov turbulence, the usual measure of 
turbulence strength is the structure constant, Cj, of the 
temperature fluctuations (Tatarskii, 1971), and the concept can 
be extended to anisotropic turbulence as well (»foo et al. , 
1980a). Since the structure constant depends upon the spectral 
index of the turbulence (which, as we have pointed out, is 
uncertain for large-scale layered atmospheric inhomogeneities), 
we will consider the general class of isotropic and layered 
media in which the resultant rms phase shift upon passage 
through the atmosphere is dominated by the largest scale 
features, with characteristic dimension Lq* The Kolmogorov 
spectrum satisfies this condition. From an obvious 
generalization of Eq. 1, we have 

< n 2>>/2 . 

The effective path through the atmosphere, L, is given by 
(2 xHH)^/ 2, and for isotropic turbulence, there will be N ■ L/Lq 
transected independent turbules of size L^. Then adopting the 
usual random-walk assumptions, we can approximate Eq. (lb) as 
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- iE ^ (IV) 

or 


a result obtained by HJW. In terms of the rms fractional 
refractivity fluctuations, <( and the normalized 
stellar intensity, ‘, we find: 

^ iso 


(2.)V* hV" A . l\//^\Al/J (19) 

Tf2 Rl/* ( ♦ / YaJ/^®° 


If the turbulence is Kolmogorov, the fractional refractivity 
fluctuations can be related to the structure constant 


(Haugstad, 1981): 

^2 1/2 a ^/3 ct (20) 

''A 2l/2(Y-l)To 

where y is the ratio of specific heats and Tq is the mean 
temperature of the background atmosphere. Then 
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!/■* W3/** hV- U5/6 
(Y-1) "^f^ 



( 21 ) 


With Eq. 13, we can relate the turbulence strength. A, 
directly to the structure constant: 

A = 2^/^ir / _H Y Cx (isotropic) (22) 

(Y-l)To Irf; 

Analogous results for layered atmospheric structure can be 
found with the modification that N, the number of Independent 
turbules crossed by a ray, is now given by H/L^: 
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or 

<*t2>i/2 ^ 2^/^ K LqS/6 fl ^ a q (24) 

(T-l) '■f^ To 

Recall that A » <*^ 2 >l /2 jqj. case ^ “ 0.5. 

Comparing Eqs. 19 and 13, we find the ratio of the rms 
phase shift for layered and isotropic turbulence with the same 
outer scale and rms refractivity to be 

(25) 

<* 2>l/2 \ H J 

t ISO 


The implication is that anisotropic turbulence should give 
enhanced scintillations compared to isotropic turbulence, for 
the usual circumstance, R/H >> 1, Physically, this is because 
there is less cancellation of phase: fewer independent 

turbules are encountered in a layered medium, and the 
contribution of each turbule is enhanced by the long slant path 
of the ray through the fluctuation. 

To obtain a feeling for the magnitudes involved, we note 
that, for Kolmogorov turbulence, 

Ct “ AT (26) 

( Ar ) ^/ ^ 

where AT is the rms temperature difference over a separation 
distance of Ar. For Uranus, we adopt * H = 70 km, R * 26000 
km, T =» 150 K, Y = 1.47, and rf = 2.45 km. Then, from Ekjs. 21 
and 29 we have 

A = 30 AT (isotropic) (27) 

(Ar)l/3 

and 


A 


I 


;S 
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A - 210 AT (layered) (28) 

(Ar)l/3 

for Ar in km. (Note that the turbulence amplitude for a 
layered atmosphere is ~7 times larger than for the isotropic 
case. ) Since the simulations show that intense scintillations 
occur for A > 1, we conclude that, in a layered atmosphere, 
temperature variations as small as 0.005 K over an altitude 
range Ar ■ 1 km can produce large spikes. 

It is unlikely that such small temperature fluctuations 
could be maintained coherently over the large horizontal 
distances (>140 km) for which individial light curve spikes 
remain recognizable (French et , 1982b). Of course, Eqs. 2^ 
and 28 assume Kolmogorov scaling of the variations. An 
alternate possibility is that the anisotropic turbulence 
spectrum is flatter than the Kolmogorov law, but Crane (1980) 
shows that in fact the vertically directed spectrum has a 
steeper slope than the Kolmogorov spectrum. Still another 
possibility is that the spikes are produced by relatively 
large-scale structure (with typical length scale > rf). We 
discuss these possibilities in greater detail in the next 
section. [Note: this transition needs to be improved]. 

III. Atmospheric Waves. Spikes, and Ray Crossing 

We have demonstrated that intense scintillations can occur 
in occultations by atmospheres with C^ • 0.01 K km"^/^. For 
comparison, a typical terrestrial value is Ct ■ 0.2 K km“^/^ in 
the lower stratosphere (Crane, 1980), and Cj * 1.5 K km”^/^ for 


Venus at 6U km altitude (Woo and Ishlmaru, The strong 

evidence for layering In the Uranus makes such comparisons 
misleading, however. It seems far more likely that the 
Uranian upper atmosphere has a combination of small-scale 
Isotropic turbulence, moderately anisotropic turbulence at 
larger scales, and underlying atmospheric structure of 
significant horizontal extent. Under these more complicated 
circumstances, what does an occultatlon curve look like? How 
well can the numerical inversion method (which assumes 
horizontal coherence over a distance L - (2wHH)^/2 (French ^ 
al . , 1978) retrieve the true atmospheric basic state in the 
presence of turbulence and horizontal Inhomogeneities? To 
explore some of these issues, we have constructed light curves 
for a variety of model atmospheres with large-scale non- 
isothermal features. 

A. Waves in the Presence of Turbulence 

If the spikes in the Uranus observations (Fig. 11) are in 
fact turbulence-induced scintillations, is it still possible 
for large-scale layered structure in the atmosphere to be 
correctly determined by numerical inversion? As a test, we 
generated, using wave optics, light curves for a model 
atmospheres with a vertical wave of constant fractional 
refractivity c (with e << 1) and vertical wavenumber u ■ 

2wn/HQ, where n is the number of vertical wavelengths per scale 
height and Hq is the mean atmospheric scale height. The 
corresponding perturbation phase component on the thin screen. 
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ORIGINAL FAGS IS 

OF POOR QUALITY 


obtalred by integrating along the ray-path through the 
atmosphf?re, is given by 

♦w(r) « sln[ ^ / 2 _tan-;^|iHQ) (29) 

•w(r) . [1 + (uHo)2]1/‘^ 

vl:ere specifies the phase of the atmospheric wave at the 
halt- light altitude, r*ro« The temperature profile 
corre‘~.;'Onciing to the sinusoidal refractivity profile is 

n by integrating the hydrostatic equation and using the 
peiiect gas law, with the result: 

H(r) = HqJ ^ + ^ 

I (1 + e sin Q) 

where the temperature is expressed in terms of scale height, 
and 

n = y(r-ro) + ♦ (31) 

In the short wavelength limit, yrio >> 1, the temperature 
profile is sinusoidal with amplitude e. In the long wavelength 
limit, yrio << 1, and the temperature wave is strongly 
suppressed, with amplitude e(yHo)^ << e. 

As a typical example, we present results for the case, n * 
3, c - 0.005. The light curve corresponding to this wave is 
shown in Fig. 12a. The stretching of the sinusoidal variations 
is clearly evident deep in the light curve. Next, we added the 
wave perturbation phase front, ♦w# to the appropriate 
perturbation phase front for strong scintillation (A “ 4.0), 
and constructed the corresponding light curve (Fig. 12b). Also 
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Shown is the light curve for turbulence alone (Fig* 12c)* 

The effects of the wave are barely perceptible by eye when 
strong turbulence is present. The results of numerical 
inversion of these light curves are shown in Pig. 13* The 
sinusoidal variation in temperature is retrieved exactly, in 
the absence of turbulence (Fig* 13a), and is still clearly 
discernible when strong scintillation is present (Fig. 13b)* 

The spurious temperature variations due to the turbulence alone 
are shown in Fig. 14c, Here, even for strong scintillations, 
these spurious variations are only a few tenths of one per 
cent. We conclude that the horizontally non-varying component 
of the thermal structure can still be retrieved, even when 
strong turbulence is present. It is intoresting that 
observations of oceanic eddies suggest that such mixtures of 
turbulence and waves do occur in fluids of geophysical interest 
(Rhines, 1H73). We note parenthetically that the rms amplitude 
of the fluctuations is almost identical for the two turbulent 
light curves in Fig. 12, so that such statistical calculations 
act as high-pass filters, and are not sensitive to low- 
amplitude long wavelength atmospheric structure. 

B. Occultations by Atmospheres With Gravity Waves 

The simultaneous presence of turbulence, gravity waves, 
and global-scale layered structure in real atmospheres 
introduces both complexity and uncertainty into the spectrum of 
refractivity variations. We have determined that the radially 
symmetric component of the atmospheric structure can be 
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retrieved In the midst of turbulence, as long as ray crossing 
is not severe, but many layered features in atmospheres are 
patchy and localized, and are not coherent over the distance L* 
( 2 wRH) 1/2 required by the inversion method. 

Occultation temperature profiles are frequently 
characterized by wavelike variations which, for the 6 Sco 
occultation by Jupiter, are consistent with gravity waves 
(French and Giernsch, ia74). The confirmation of this 
interpretation would be significant, since propagating inertia 
gravity waves can be important to the upper atmospheric energy 
balance. To explore this issue, we have constructed model 
light curves for an isothermal atmosphere with a realistic 
field of gravity waves, using geometric optics. We assume a 
vertical wavelength Xy ■ H, and horizontal wavelengths X^ of 
20, 15 and 10 times the scale height, H. The radius of the 
planet is chosen to be R - 1000 H. The ratio, X^/L, is a 
measure of the horizontal coherence of the atmosphere. For 
X^/L < 1, Che assumptions of the inversion method are no longer 
satisfied. 

In the absence of damping, temperature perturbations grow 
exponentially with height, for gravity waves, but in real 
atmospheres such growth is checked by internal instabilities, 
or by trapping and reflection of the waves. We model the 
situation where the eddy diffusivity is just sufficient to 
suppress the exponential growth, so that the fractional 
amplitude, c, of the wave is constant with height; we take c ■ 
0.05. The refractivity field is definied, after French and 
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Glerasch (1974), as 

v(r, 2 ) - + c sin lii(r-ro) kz + t)J^ (32) 

where m and k are the vertical and horizontal wavenumbers, 
respectively. 

The resulting light curves are displayed In Pig. 14. The 
top figures corresponds to a radially symmetric atmosphere. As 
Is decreased, the horizontal structure changes significantly 
over L, the effective path length through the atmosphere, and 
the broad spikes in the corresponding light curves are 
suppressed and, ultimately, become Imperceptible. The 
corresponding temperature profiles, obtained by inversion of 
the light curves, are shown in Fig. 16. For X^/L < 0.25, the 
vertical structure of the atmosphere can no longer be 
determined accurately by numerical Inversion because of 
horizontal averaging effects. 

C. Severe Hay Crossing 

The preceding simulations demonstrate that horizontal 
inhomogeneities produced by either turbulence or by larger 
scale waves have little effect on the inferred temperature 
profile, and that the uuasl-global component of the atmospheric 
structure can be accurately retrieved even when turbulence and 
short-wavelength waves are strong. This is not true in 
general: an implicit assumption of the inversion technique Is 

that ray crossing is not severe, since the method requires that 
the light at a point on the light curve corresponds to a unique 
ray. In the wave model for spikes, it can be shown that the 
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condition for eventual ray crossing, as 4 -*^0, is: 


e(wH)V2 > 1 

(33) 

e(2tn)3/2 > 1 

(34) 


Thus, for a given vertical wavelength, there is a critical wave 
amplitude, e, above which ray crossing will occur. Fig. 16 
displays the regimes of no ray crossing and of eventual ray 
crossing. 

What is the effect of ray crossing on the inversion 
technique? In any inversion, the light curve is divided into 
bins and the average flux in computed for each bin. If the ray 
crossing is sufficiently mild to be confined to a single bin, 
then the inversion will not be affected, since the mean flux 
does not depend on how the rays have distributed themselves in 
any given bin. We find that the minimum altitude interval that 
can be resolved. Ah, is given by 

^ (35) 

H 

or 

Ah > c(2wn) (36) 

H 

In order to resolve a wave of n wavelengths per scale height, 
we require Ah < V^nH; thus, the maximum amplitude wave that 
can be resolved when ray crossing is present is given by 

c < l/(8«) l/2nV2 (37) 

Two model light curves with ray crossing are shown in Fig. 
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17. In the first case, a fairly large vertical wavelength 
(i.. » H/2, or n “ 2) of amplitude e ® 0. 1 is assumed. 

According to Eq, 37, the maximum amplitude wave that can be 
resolved is e = 0.07, so that ray crossing is severe. In the 
second case, Xy * H/5, n = 5, and e = 0.025. The critical 
amplitude from Eq. 37 is c = 0.018, and again, ray crossing is 
severe. In each case, the refractivity perturbations produce 
intense, sharp spikes in the light curve. Inversion of the 
light curves yields the derived temperature profiles shown in 
Fig. 17. At first, the correct atmospheric structure is 
derived. As the inversion proceeds and ray crossing becomes 
severe, the computed profile shows a strongly damped wave, 
whereas in fact the true profile has a constant amplitude wave. 
The result is that the inversion essentially "ignores" the 
severely crossed rays, and, importantly, introduces only very 
minor spurious features into the derived profiles. When the 
ray crossing is not severe, the inversion technique introduces 
insignificant erro. into the retrieved profiles, a result also 
obtained by Wallace (1975). 

All of the foregoing simulations tend to confirm the 
following assertion: Quasi-global features of atmopsheric 

structurp are accurately determined by numerical inversion ; 
horizontally inhomogeneous structures (including turbulence) 
with coherence scale L (2irRH) ^ ^ are filtered out and have 
little effect on the derived temperature profiles. Naturally, 
the influence of photon noise and terrestrial scintillation 
must also be considered when interpreting the observations. 
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D. The Wave Optics of Ray Crossing 

The light curves In Fig. 17 were computed using geometric 
optics, which is not legitimate for determining the detailed 
structure of well-focussed or over-focussed spikes. We have, 
therefore, computed light curves of a sequence of Increasingly 
sharp spikes to determine the effects of diffraction in the 
ray-crossing regime. In Fig. 19, a sequence of spikes is 
displayed. The upper figures were computed with wave optics, 
and for the lower figures we employed geometric optics. In 
each case, an atmospheric wave with vertical wavelength Xy “ 

H/4 was used. Each figure shows the same region of the 
resultant light curve; the only differences are in the 
amplitude of the wave, e, and in the optics assumed. (Note, 
however, that the vertical scale varies from plot to plot.) 

For the first case (e “ 0.01), wave optics and geometric optics 
give essentially identical results. For c = 0.015, diffraction 
fringes become apparent for this very sharp spike, and for the 
ray-crossed spikes in the two final cases, there are strong 
interference fringes in the wave optical results. 

Although the detailed structure of the spikes clearly 
depends on the assumed atmospheric structure, a general result 
is that a single layered feature can be responsible for many 
detailed variations in the light curve. This suggests that 
multiple-peaked "blobs" along the tail of an occultation light 
curve may represent diffraction of light passing through a 
single coherent layer, rather than through multiple layers. 

(An example of such features is shown in Fig. 2c of French 
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et al^. (1982b), for the 15 August 1980 occultation by Uranus.) 
IV. Conclusions 

Many of the problems of interpretation of stellar 
occultation observations stem from the difficulty of defining 
in a general way the effects that realistic atmospheric 
structure have on the resultant light curves. The theory of 
turbulent scintillations in occultations has been well 
developed, for conditions of interest, only for weak 
scintillations (HJW), whereas the assumption of perfect 
layering in the atmosphere has permitted temperature profiles 
to be derived whose validity is uncertain. We have formulated 
quite general techniques for producing model light curves for a 
wide variety of realistic atmospheres, not subject to the above 
restrictions. We generally confirm the weak-scintillation 
results of HJW, but we find that there are important 
qualitative differences for strong scintillations. These 
include saturation of the modulation index at unity for strong, 
one-dimensional, Kolmogorov turbulence, and independence of the 
FWHM of the ACF of light curve residuals upon the mean 
intensity, in strong contradiction to weak scintillation 
theory. Comparison of recent observations with the turbulence 
simulations reveals that the modulation index of the 
observations saturates at unity, suggestive of strong 
scintillation, but that the FWHM of the ACF is significantly 
wider than the turbulence simulations. The critical test of 
the importance of diffraction effects will be to observe an 
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occultatlon cf a small star at two widely separated 
frequencies. 

Although It remains an open question whether or no«. the 
fine-scale structure observed during occultatlons is turbulent 
scintillation, we have demonstrated In a variety of ways that 
the numerical inversion of light curves is only weakly 
influenced by horizontally Inhomogeneous atmosperic structure, 
and that the underlying large-scale coherent structure can be 
retrieved. This suggests the possibility that occultatlons 
will allow us to determine both the quasi-global atmospheric 
structure and the statistical characteristics of small-scale 
refractivity variations. A large number of high-quality 
observations of occultatlons exist, and a most productive 
approach would be to construct model turbulent light curves for 
the specific parameters (H, rf , and d») appropriate to each 
event. By comparing the statistical properties of the 
observ:it ions with the simulations, it should be possible to 
restrict considerably the possible turbulence characteristics 
in these planetary upper atmospheres. 
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APPENDIX 1 

Turbulence Amplitude 
Autocorrelation Function 
Generalized Correlation Function 
Structure Constant of Turbulence 
Cross-Correlation Function 
Earth-Planet distance 
Projected Stellar Diameter (km) 

Non-dimensional Stellar Diameter (d*/rf) 

Point Source Wave Amplitude at the Earth 

Point Source Wave Amplitude at the Reference Screen 

Fast Fourier Transform 

Atmospheric Scale Height 

Non-Dimensional Scale Height (H/rf) 

Scale Height of Mean Atmosphere 
Point Source Wave Intensity at the Earth 
Finite Source Wave Intensity at the Earth 
Isotropic 

Horizontal Wavenumber of Inertia-Gravity Wave 
Effective Path Length Through the Atmosphere 
[(2 itRH)1/2] 

Outer Scale of Turbulence 

Modulation Index of Light Curve Fluctuations 
Number of Turbules Transected by a Ray 
Number of Vertical Wavelengths per Scale Height of an 
Atmospheric Wave 

Wavenumber Power Spectrum of Phase on the Reference 
Screen 

Wavenumber Power Spectrum of Atmospheric Refractivity 
Fluctuations 

Vertical Wavenumber of Turbulent Refractivity 

Fluctuations 

Radius of Planet 

Distance from Center of Planet 

Fresnel Scale [(XD)^/2] 

Distance from Center of Planet of Half-light Ray 
Brightness Distribution of the Star 
Half-width of Interval Over Which Correlation is 
Computed 

Mean Temperature of Background Atmosphere 
Vertical Coordinate 

Horizontal Coordinate Perpendicular to Earth-Planet 
Line 

Horizontal Coordinate Parallel to Earth-Planet j-j’-is 
tvaiiu of Specific Heats 
Altitude Interval in Atmosphere 
Interval in Altitude 

RMS Temperature Difference Over a Separation Distance 
of Ar 
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Deviation of Normalized Light Curve Intensity from 
Mean Value, 

e Fractional Amplitude of Atmospheric Wave in 

Ref ractivity 

Ox Subtended Angie in x Direction as Viewed from Earth 

8y Subtended Angle in y Direction as Viewed from Earth 

6* Subtended Angle of Star as Viewed from Earth 

X Wavelength of Light 

Xji Horizontal Wavelength of Inertia-Gravity Wave 

Xy Vertical Wavelength of Inertia-Gravity Wave 

ji Vertical Wavenumber of Inertia-Gravity Wave 

V Atmospheric Refractivity 

VA Mean Atmospheric Refractivity 

vt Turbulent Atmospheric Refractivity 

Vq Atmospheric Refractivity at Altitude of Half-light 

Ray 

p Axial Ratio of Anisotropic Turbules (Horizontal/ 

Vertical ) 

T Lag of Correlation Function 

♦ Phase Shift of Plane Wave Produced by Passage Through 

the Atmosphere 

♦a Contribution of Mean Atmosphere to ♦ 

Contribution of Turbulence to ♦ 

Contribution of Atmospheric Waves to ♦ 

4 Mean Intensity of Light Curve, In Units of Intensity 

of the Unocculted Star 

\(i Phase of Atmospheric Wave at Half-light Altitude, rQ 

Q Phase of Atmospheric Wave at Altitude, r 
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FIGURE CAPTIONS 

Fig. 1. Occultation Geometry* The distortions in the wave 

front of a plane wave, due to passage through a turbulent 
atmosphere, can be represented by a thin phase-changing 
screen at distance D from the earth's path. The light 
curve as viewed from the earth along its path is computed 
in overlaping sections (for example, A and B) by 
summing the phases of the appropriate section of the 
upstream phase screen. Edge effects are avoided by 
apodizing. 

Fig. 2. Occultation light curves for model isothermal 

atmospheres with anisotropic Kolmogorov turbulence. A 
measure of the turbulence strength is given by A, the rms 
turbulence-induced phase change (in radians) at the level 
in the atmosphere corresponding to the half intensity ra’' 
For A ^ 1, the resultant scintillations are strong. 

Fig. 3. Turbulent occultations of a finite source. In the 
simulations at left, a weak-scintillatlon (A ■ 0.5) 
point-source light curve is convolved with succesively 
larger sources whose diameter at the occulted planet is d* 
(in units of rf). On the right, a strong scintillation 
(A » 4.0) case is shown. The scintillations are strongly 
damped for d* ^2 rf. 

Fig. 4. Statistical properties of turbulent scintillations. 

The rms amplitude of the fluctuations from the mean of the 
point-source light curves in Fig. 2 is plotted as a 
function of distance along the earth's path. The 
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triangles mark the measured values and the solid line Is 
from HJW's weak scintillation theory. Also shown is an 
Isothermal light curve. Note that the scales on 
successive plots are Increased by a factor of two. As the 
turbulence strength Increases beyond A • 2, the 
scintillations saturate, a phenomenon not predicted by 
weak scintillation theory. At right, the modulation index 
Is shown for the same circumstances, along with the weak 
sclntllla.lon model of HJW. As the turbulence becomes 
very strcig, the modulation Index saturates to unity. 

Fig. 5. Autocorrelation function of point source light curve 
residuals for weak (A - 0.5) and strong (A ■ 4.0) 
turbulence. For weak scintillation, the ACF broadens as 
the stellar flux, decreases, but for strong scintill- 
ation, the ACF remains essentially Independent of 4 . 

Pig. 6 . Point source autocorrelation FWHM for residuals of 
turbulent light curves. For weak turbulence, the weak- 
sclntlllatlon model, shown by a bold solid line. Is 
obeyed: the FWHM Increases as 4 ”^/ 2, For turbulence 

amplitudes A > 2, the FWHM of the ACF remains small at 
all 4 . 

Fig. 7. Two-frequency cross-correlation of turbulent 

point-source light curve residuals for a frequency ratio 
of 2 . The prediction of HJW's weak scintillation theory 
Is shown dashed. As the turbulence amplitude Is 
increased, the maximum cross-correlation diminishes. 

Fig. 8 . Theoretical power spectrum of light curve residuals 
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for weak anisotropic turbulence, calculated from HJW's 
model. See text for details. 

Fig. Computed power spectra of light curve residuals for 
model turbulent occultations of a finite source with 
projected diameter d* (in units of Tf) and turbulence 
amplitude A ■ 0.5, corresponding to weak scintillation. 

Pig. 10. Computed power spectra of light curve residuals for 
model turbulent occultations of a finite source with 
projected diameter d* (in units of rf) and turbulence 
amplitude A > 4.0, corresponding to strong scintillation. 

Fig. 11. Occultations observations from CTIO of the 8/15/80 
Uranus occultatlon of KM 12. The immersion light curve 
contains many sharp spikes. The rms amplitude and 
modulation index of the fluctuations are plotted as 
triangles, below. Also shown is the model isothermal 
light curve scaled appropriately for the occultatlon. 

Fig. 12. Model occultatlon light curves for radially symmetric 
atmospheres containing waves and/or strong turbulence, 
computed using wave optics. Notice that, in the presence 
of turbulence, the effects of the wave are barely 
perceptible to the eye. 

Fig. 13. Temperature profiles obtained by numerical inversion 
of the light curves in Fig. 12. The wave is still clearly 
resolved, in the presence of strong turbulence. 

Fig. 14. Model occultatlon light curves for atmospheres 

containing inertia-gravity waves. See text for details. 

Fig. 15. Temperature profiles obtained by numerical inversion 
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Of the light curves in Fig. 14. As the horizontal 
wavelength of the gravity waves is decreased, the 
inversion technique is less and less able to determine 
the true atmospheric basic state. 

Fig. Iti. Ray-crossing regimes. 

Fig. 17. Model light curves, computed from geometric optics, 
for two cases of severe ray crossing. Intense spikes are 
caused by the focussing of starlight by the atmospheric 
waves. 

Fig. Id. Temperature profiles obtained by numerical inversion 
of the light curves in Fig. 17. Deep into the 
occultation, severe ray crossing prevents the true 
atmospheric structure from being determined. 

Pig. 19. Comparison of light curve spikes computed from wave 
optics (upper) and geometric optics (lower) for 
atmospheric waves whose fractional amplitude, c, is 
increased from 0.01 (left figures) to 0,025 (right 
figures). When diffraction effects are taken into 
account, the light curve contains many fringes. Note the 
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TURBULENCE MODELS AND THE ATMOSPHERE OF URANUS 
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APPENDIX F 


STRUCTURE OF THE UPPER ATMOSPHERE OF URANUS 


DEPARTMENT OP EARTH AND PLANETARY SCIENCES 


MASSACHUSETTS INSTITUTE OF TECHNOLOGY 

CAMBRIDGE, MASSACHUSETTS 02139 


Abstract: The 1980 August 15 Occultation of KM 12 by 

Uranus as Observed Frotii CTIO: Atmospheric 

Results 

Richard G. French, J.L. Elliot, et al. 


Observations from CTIO of the occultation of KM 12 
by Uranus on 1980 August 15 have been reduced to obtain 
temperature profiles of the Uranun upper atmosphere for 
number densities between 10'* and 10“ cm“*. The profiles 
are characterized by wavelike variations of i5 K and 
vertical wavelength of two scale heights. The deduced 
variations are remarkably similar to those obtained from 
occultation observations on 10 March 1977 and 10 June 1979, 
suggesting the existence of a long-lived guasi-global 
thezrmal layer in the Uranus upper atmosphere. Radiative 
time constants are long at these altitudes, but smoothing 
of temperature fluctuations by eddy mixing and molecular 
diffusion is rapid. Consequently, it is unlikely that 
direct deposition of solar radiation (as in a Chapman 
layer) is responsible for the observed structure. The 
maintenance of this non- isothermal layer requires a signif- 
icant source of energy which is important to the upper 
atmosphere heat balance. 


